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ABSTRACT 

The  measurable  quantities  associated  with  7-ray  and  neutron  observations  of  solar  flares  are  nuclear-deexcitation 
line  shapes,  shifts,  fluences,  and  time  histories;  neutron  capture  and  annihilation  line  fluences  and  time  histories;  and 
energy-dependent  escaping  neutron  fluence  and  time  history.  A  comprehensive  understanding  of  these  quantities 
requires  a  model  for  ion  acceleration,  transport,  and  interaction.  In  this  paper  we  address  transport  and  interaction 
using  a  magnetic  loop  model  that  includes  energy  losses  due  to  Coulomb  collisions,  removal  by  nuclear  reactions, 
magnetic  mirroring  in  the  convergent  flux  tube,  and  MHD  pitch-angle  scattering  in  the  corona.  The  accelerated  ions 
are  assumed  to  have  a  given  kinetic  energy  spectrum  and  composition.  Each  measurable  quantity  depends  to  varying 
degree  on  the  parameters  of  the  loop  model  and  of  the  accelerated  ions.  We  explore  these  dependences  in  detail  and 
construct  a  self-consistent  approach  to  the  analysis  of  high-energy  flare  data  that  provides  an  optimum  set  of  param¬ 
eters  with  meaningful  uncertainties.  To  illustrate  this  approach,  the  calculations  are  applied  in  a  comprehensive 
analysis  of  the  7-ray  and  neutron  observations  of  the  1991  June  4  solar  flare  obtained  with  OSSE  on  CGRO.  We  find 
that  the  loop  model  can  account  for  these  observations  with  physically  reasonable  values  for  the  parameters.  In  addition, 
our  analysis  of  the  neutron  data  shows  that  the  accelerated  ion  spectrum  for  this  flare  was  not  an  unbroken  power  law 
but  had  to  steepen  sharply  above  ~125  MeV  nucleon-1.  The  paper  also  provides  yields  and  yield  ratios  calculated 
with  assumed  abundances  and  spectral  forms  currently  considered  appropriate  for  solar  flares.  They  can  be  used  by  other 
researchers  analyzing  high-energy  solar  flare  data. 

Subject  headings:  acceleration  of  particles  —  nuclear  reactions,  nucleosynthesis,  abundances  —  Sun:  flares  — 
Sun:  X-rays,  gamma  rays 


1.  INTRODUCTION 

The  interactions  of  ions  accelerated  in  solar  flare  magnetic  loops 
produce  excited  and  radioactive  nuclei,  neutrons,  and  7r-mesons. 
All  ofthese  products  subsequently  produce  observable  7-ray  emis¬ 
sion  via  secondary  processes  and  the  neutrons  that  escape  may  be 
observed  directly  in  space  and  on  Earth  and  indirectly  via  their 
decay  protons.  The  measurable  quantities  associated  with  these 
emissions  (e.g.,  7-ray  line  intensities,  shapes,  and  time  histories 
and  escaping  neutron  spectra  and  time  histories)  convey  informa¬ 
tion  about  the  accelerated  ions  responsible  for  their  production. 
This  information  is  described  by  acceleration  parameters,  such 
as  the  number  and  spectral  shape  of  the  accelerated  protons 
and  a-particles  (and  thus  their  energy  content)  and  the  relative 
abundances  of  the  accelerated  ions,  and  it  can  be  obtained  from 
measured  deexcitation  line  intensities  and  ratios.  This  aspect  of 
high-energy  flare  emission  has  been  well-explored  in  past  re¬ 
search  (e.g.,  Ramaty  et  al.  1995). 

However,  the  measurable  quantities  can  also  depend  strongly 
on  the  physical  conditions  of  the  magnetic  flare  loop  in  which  the 
particles  transport  from  the  acceleration  site  to  the  interaction 


region.  These  physical  conditions  are  described  by  physical  pa¬ 
rameters,  such  as  loop  length,  the  level  of  pitch-angle  scattering 
(PAS)  within  the  loop,  the  degree  of  convergence  of  the  loop 
magnetic  field  at  the  footpoints,  the  location  of  the  flare  on  the 
Sun,  and  the  density  and  temperature  height-profile  and  elemen¬ 
tal  composition  of  the  flaring  atmosphere.  Analyses  of  high- 
energy  solar-flare  emissions  therefore  offer  the  potential  to  learn 
much  more  about  the  flare  process  than  just  particle  acceleration, 
such  as  the  structure  and  evolution  of  the  flare  environment. 

But  because  each  of  the  measurable  quantities  can  depend  on 
some  or  all  of  the  same  acceleration  or  physical  parameters,  de¬ 
riving  well-constrained  values  for  those  parameters  is  challeng¬ 
ing.  Previous  analyses  often  have  not  taken  these  complicated 
interdependences  into  account  in  a  systematic  way,  relying  on  ad 
hoc  assumptions  that  resulted  in  unrealistically  small  uncertain¬ 
ties  for  their  limited  number  of  parameters.  The  goal  of  this  paper 
is  to  systematically  address  these  interdependences  by  identify¬ 
ing  all  relevant  parameters  and  showing  how  the  nature  of  the 
nuclear  interactions  depend  on  them  and  how  the  various  mea¬ 
surable  quantities  are  subsequently  affected.  This  understanding 
will  enable  us  to  construct  a  comprehensive  and  self-consistent 
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approach  to  the  analysis  of  flare  data  that  will  provide  optimum 
parameter  values  (and  meaningful  uncertainties)  that  predict  ob¬ 
servable  quantities  that  are  consistent  with  all  of  the  measured 
quantities.  Success  of  this  approach  requires  that  the  available 
flare  measurements  cover  a  sufficiently  broad  range  of  the  mea¬ 
surable  quantities.  The  included  figures  describing  the  dependence 
of  the  measurable  quantities  on  the  acceleration  and  physical  pa¬ 
rameters  can  be  used  by  other  researchers  analyzing  high-energy 
solar  flare  data. 

In  §  2  we  describe  in  detail  how  the  acceleration  and  physical 
parameters  affect  the  nuclear  interactions  responsible  for  the  high- 
energy  flare  emissions  and  how  the  measurable  quantities  asso¬ 
ciated  with  these  emissions  are  subsequently  affected.  We  also 
explore  the  ranges  of  accelerated  ion  energies  that  are  responsible 
for  the  various  emissions.  With  this  understanding,  we  develop  in 
§  3  a  data  analysis  algorithm  that  provides  the  approach  and  anal¬ 
ysis  steps  necessary  for  a  comprehensive  analysis  of  solar  flare 
7-ray  line  and  neutron  data.  As  a  demonstration  of  the  procedure, 
we  use  the  measurements  of  the  1991  June  4  flare  to  constrain  all 
of  the  relevant  parameters.  The  1991  June  4  flare  was  observed 
by  the  Oriented  Scintillation  Spectrometer  Experiment  (OSSE) 
on  the  Compton  Gamma  Ray  Observatory  ( CGRO )  and  by  a  num¬ 
ber  of  other  instruments.  The  flare  is  one  of  the  best-measured 
high-energy  flares  to  date,  and  the  resulting  data  set  is  unique 
in  that  it  provides  a  remarkably  wide  range  of  time-dependent 
measurements  with  good  statistical  significance  of  hard  X-rays, 
7-ray  lines  and  continuum,  and  neutrons  as  well  as  data  at  other 
wavelengths. 

2.  GAMMA-RAY  AND  NEUTRON  PRODUCTION 
BY  ACCELERATED  IONS  IN  SOLAR  FLARES 

Accelerated-ion  interactions  with  the  ambient  solar  medium 
produce  excited  nuclei  that  promptly  deexcite,  yielding  observ¬ 
able  7-ray  lines.  The  lines  resulting  from  accelerated  protons  and 
o-particles  interacting  with  ambient  nuclei  heavier  than  He  are 
“narrow”  (fractional  FWHM  ~2%),  due  to  the  relatively  low  re¬ 
coil  velocity  of  the  heavy  nucleus.  Accelerated  a-particles  and 
heavier  ions  interacting  with  ambient  H  and  He  produce  “broad” 
lines  (fractional  FWHM  ^20%)  due  to  the  relatively  high  recoil 
velocity  of  the  nucleus.  Recently,  Kuzhevskij  et  al.  (2005a)  con¬ 
sidered  deexcitation  line  yields  from  heavy-heavy  nucleus  inter¬ 
actions,  but  we  have  found  that  for  the  narrow  1.634  MeV  20Ne 
and  1.369  MeV  24 Mg  lines  produced  by  interactions  of  12C  and 
160,  the  yields  are  negligible  due  to  both  the  lower  abundances  of 
both  species  and  the  higher  energy  loss  of  heavy  particles.  Ac¬ 
celerated  ion  interactions  also  create  radioactive  nuclei,  which 
produce  delayed  7-ray  lines  that  can  be  observed  when  the  prompt 
emission  is  negligible  (Tatischeff  et  al.  2006). 

Ion  interactions  also  produce  neutrons  that  can  be  observed 
both  directly  and  indirectly.  Neutrons  that  escape  from  the  Sun 
and  survive  the  transit  to  Earth  can  be  directly  detected  with  in¬ 
struments  in  orbit  and,  for  events  with  sufficient  high-energy  neu¬ 
trons  (>200  MeV),  by  neutron  monitors  on  the  ground.  Those 
escaping  neutrons  that  decay  can  be  indirectly  detected  via  their 
decay  protons  in  space  (e.g.,  Evenson  et  al.  1983).  Neutrons  that 
remain  at  the  Sun  may  be  indirectly  detected  because  some  of 
them  can  be  captured  on  ambient  hydrogen  in  the  photosphere, 
producing  deuterium  and  a  strong  7-ray  line  at  2.223  MeV. 

Our  understanding  of  7-ray  and  neutron  production  has  evolved 
for  more  than  35  years.  After  the  seminal  work  by  Lingenfelter  & 
Ramaty  (1967)  a  comprehensive  treatment  of  deexcitation-line 
production  was  given  by  Ramaty  et  al.  (1979).  The  resulting  7-ray 
line  production  code  has  recently  been  updated  (Kozlovsky  et  al. 
2002)  to  incorporate  new  cross  section  data.  This  code  calculates  a 


complete  7-ray  line  spectmm  along  any  observing  direction  from 
the  flare  site  for  any  assumed  accelerated-ion  spectral  shape,  accel¬ 
erated  and  ambient  composition,  and  interacting-particle  angular 
distribution.  We  will  use  this  deexcitation-line  production  code  in 
this  paper.1 

Hua  et  al.  (1989)  developed  a  magnetic  loop  transport  and  in¬ 
teraction  model  (described  more  fully  below)  in  which  a  Monte 
Carlo  simulation  follows  an  individual  accelerated  ion  through  a 
flare  loop  until  it  either  undergoes  a  nuclear  reaction  or  its  energy 
falls  below  nuclear  reaction  thresholds.  The  process  is  repeated 
until  adequate  statistical  significance  is  achieved.  The  most  im¬ 
portant  aspects  of  ion  transport  are  included  in  the  model:  energy 
losses  due  to  Coulomb  collisions,  removal  by  nuclear  reactions, 
magnetic  mirroring  in  the  convergent  flux  tube,  and  MHD  pitch 
angle  scattering.  They  applied  this  model  to  deexcitation  line  re¬ 
actions  and  calculated  the  depth  distribution  and  time  history 
(for  instantaneous  release  of  the  energetic  ions)  for  production  of 
several  nuclear  deexcitation  lines.  We  will  use  this  deexcitation¬ 
line  transport  code  in  this  paper. 

The  theory  of  neutron  production,  escape,  and  transport  to 
Earth  was  first  treated  by  Lingenfelter  et  al.  ( 1 965a,  1 965b)  and 
was  further  developed  by  Hua  &  Lingenfelter  (1987a)  using 
Monte  Carlo  transport  simulations  to  follow  the  neutrons  until 
they  decay,  are  captured  on  H  in  the  solar  atmosphere  to  produce 
the  2.223  MeV  neutron-capture  line,  or  escape  from  the  Sun.  Us¬ 
ing  the  magnetic  loop  model  referred  to  above,  Hua  et  al.  (2002) 
calculated  angle-dependent  neutron  spectra  using  updated  neu¬ 
tron  production  cross  sections  and  new  kinematics.  The  code 
calculates  escaping  neutron  spectra  and  neutron-capture  line 
intensities  and  time  histories  along  any  observing  direction  from 
the  flare  site.  We  will  use  this  neutron /neutron-capture  line 
transport/production  code  in  this  paper.  The  code  is  available 
online  (see  Hua  et  al.  2002). 

Ion  interactions  also  produce  radioactive  nuclei  and,  if  there 
are  sufficient  very  high-energy  accelerated  ions,  positively  charged 
pions,  both  of  which  decay  to  produce  positrons,  which  then  an¬ 
nihilate  to  produce  the  strong  0.511  MeV  annihilation  line.  Pro¬ 
duction  of  positron-emitting  nuclei  was  discussed  by  Kozlovsky 
et  al.  ( 1 987, 2004)  and  Tatischeff  et  al.  (2006).  Pion  production  was 
discussed  by  Murphy  et  al.  (1987).  Positron  annihilation  in  the  so¬ 
lar  atmosphere  was  recently  addressed  by  Murphy  et  al.  (2005). 
Unfortunately,  the  intense  flux  of  X  rays  from  the  June  4  flare 
restricted  the  usable  data  to  those  from  the  detector  pointed  away 
from  the  Sun  (see  §  3.1).  Because  of  the  resulting  strong  attenu¬ 
ation  of  <  1  MeV  photons  by  the  detector  shields,  measurements  of 
the  annihilation  line  during  the  peak  of  the  flare  are  not  available 
and  positron  production  will  not  be  discussed  in  this  paper. 

Ion  interactions  can  also  produce  pions.  Neutral  pions  decay 
into  two  7-rays  with  energies  of  67.6  MeV  in  the  pion  rest  frame. 
The  energies  of  these  photons  are  then  Doppler-shifted  by  the  pion 
motion.  The  spectrum  is  a  broad  peak  centered  near  70  MeV  with 
an  FWHM  of  ^  1 00  MeV  ( Lingenfelter  &  Ramaty  1 967 ;  Crannell 
etal.  1979;  Murphy  etal.  1987).  The  positrons  and  electrons  from 
charged  pion  decay  produce  continuum  bremsstrahlung,  most  eas¬ 
ily  observed  at  >1 0  MeV.  We  will  see  that  the  accelerated  particle 
spectmm  responsible  for  7-ray  emission  in  the  June  4  flare  was 
relatively  steep  and  did  not  produce  significant  numbers  of  pions. 
Pion  production  will  also  not  be  discussed  in  this  paper. 

In  §  2. 1  we  describe  the  loop  transport  and  interaction  model 
and  its  parameters  relating  to  the  physical  conditions  of  the  mag¬ 
netic  loop.  In  §  2.2  we  summarize  the  acceleration  and  physical 

1  A  simplified  version  of  the  code  is  available  at  http://lheawww.gsfc.nasa 
.  go  v/users/ramaty/  ViewPubs/ ramaty.  html . 
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parameters  addressed  by  our  calculations  of  deexcitation  7-ray, 
neutron,  and  neutron-capture  line  production.  In  §  2.3  we  discuss 
how  these  parameters  affect  the  nuclear  interactions.  (Some  of 
the  material  concerning  the  interaction  time  history  and  depth 
distribution  was  first  discussed  by  Hua  et  al.  1989  but  for  a  Bessel 
function  accelerated-particle  spectrum  rather  than  the  power-law 
spectrum  presented  here.  Some  of  the  material  concerning  neutron 
production  was  discussed  by  Hua  et  al.  2002,  but  they  did  not  in¬ 
clude  a  discussion  of  the  neutron  capture  line.)  In  §  2.4  we  describe 
how  the  measurable  quantities  associated  with  the  products  of 
these  interactions  (deexcitation  lines,  neutrons,  and  the  neutron- 
capture  line)  are  subsequently  affected  and  the  ranges  of  ion  ener¬ 
gies  that  are  most  effective  in  producing  each.  The  understanding 
gained  from  the  single-loop  study  can  be  easily  extended  to  ac¬ 
count  for  multiple-loop  structures. 

2.1.  Magnetic  Loop  Transport  and  Interaction  Model 

The  magnetic  loop  will  form  the  basic  geometrical  structure  for 
the  calculations  performed  here.  Imaging  observations  by  Yohkoh 
and  RHESSI  indicate  that  X-  and  7-ray  flare  emission  occurs  at 
footpoints,  suggesting  the  presence  of  a  loop.  Observations  at 
radio,  optical,  and  soft  X-rays  clearly  reveal  the  loop  structure.  We 
will  perfonn  our  calculations  for  a  single  loop  since  the  funda¬ 
mental  processes  probably  occur  within  one  loop.  Complex  flare 
processes  may  be  considered  as  resulting  from  combinations  of 
loops;  e.g.,  recent  imaging  results  from  RHESSI  show  that  the 
electron  bremsstrahlung  and  neutron  capture  line  emissions  are 
physically  separated,  suggesting  at  least  two  loops  (Hurford  et  al. 
2003). 

We  will  use  the  loop  formalism  developed  by  Hua  et  al.  ( 1 989). 
Their  transport  and  interaction  model  consists  of  a  semicircular 
coronal  portion  of  length  L  and  two  straight  portions  extending 
vertically  from  the  transition  region  through  the  chromosphere 
into  the  photosphere.  Recent  7-ray  line  shift  measurements  by 
RHESSI  (Smith  et  al.  2003)  suggest  that  in  some  flares  the  loop 
may  be  inclined  to  the  solar  surface.  Here  we  will  only  perform 
calculations  for  vertical  loops.  Qualitative  implications  for  inclined 
loops  can  be  estimated  from  the  results  for  a  vertical  loop. 

Below  the  transition  region,  the  magnetic  field  strength  is  as¬ 
sumed  proportional  to  a  power  6  of  the  pressure  (Zweibel  & 
Haber  1983).  For  coronal  and  photospheric  pressures  of  0.2  and 
105  dyn  cm-2  and  relevant  magnetic  field  strengths  of  100  and 
1600  G,  respectively,  b  ~  0.2.  Such  a  converging  magnetic  field 
can  result  in  mirroring  of  accelerated  particles.  Pitch  angle  scat¬ 
tering  is  characterized  by  A,  the  mean  free  path  required  for  an 
arbitrary  initial  angular  distribution  to  relax  to  an  isotropic  distri¬ 
bution.  The  dependence  of  A  on  particle  energy  is  expected  to  be 
weak  (see  discussion  by  Hua  et  al.  1989)  and  is  assumed  here  to 
be  independent.  The  level  of  PAS  is  characterized  by  2,  the  ratio 
of  A  to  the  loop  half-length  Lc  (=L/ 2).  Magnetic  convergence 
and  PAS  directly  affect  the  angular  distribution  of  the  accelerated 
particles  when  they  interact  with  the  ambient  medium.  As  we 
will  see,  this  angular  distribution  will  have  a  strong  effect  on  the 
subsequent  observable  emissions.  In  the  past,  researchers  have 
often  assumed  specific  angular  distributions  (such  as  isotropic, 
downward  beam,  fan  beam,  etc.).  With  the  present  calculations, 
the  distributions  will  develop  naturally  from  the  physical  param¬ 
eters  of  the  loop  model  and  cover  a  wider  range  of  distributions. 

Because  recent  results  based  on  measurements  made  with 
RHESSI  (Share  et  al.  2004;  Murphy  et  al.  2005;  Aschwanden 
et  al.  2002)  and  the  Solar  Maximum  Mission,  Yohkoh,  and  Granat 
Phebus  (Kuzhevskij  et  al.  2005b)  have  suggested  that  the  flaring 
atmosphere  within  the  loop  may  be  considerably  different  from 


TABLE  1 
Parameters 


Type  Parameter 


Acceleration .  Power-law  index,  s 

Acceleration  release  time  history,  axon(t) 
Accelerated-particle  composition 

Physical .  Level  of  PAS,  A 

Magnetic  convergence,  8 
Loop  half-length,  Lc 
Ambient  composition 

Atmospheric  density  and  temperature  model,  n(h)  and  T(h) 
Flare  heliocentric  angle,  0Obs 


that  of  the  quiet  Sun,  several  height-density  profiles  («(/?))  for  the 
solar  atmosphere  will  be  considered. 

In  the  Monte  Carlo  simulation,  accelerated  ions  are  released 
instantaneously  and  isotropically  at  the  top  of  the  loop.  Each  ion 
is  followed  until  it  interacts  or  thermalizes,  typically  near  the 
loop  footpoints  in  the  chromosphere  or  upper  photosphere  where 
the  density  is  sufficient  for  efficient  nuclear  interactions.  When 
an  interaction  occurs,  the  time  after  the  initial  release  of  the  ener¬ 
getic  ion,  the  location  along  the  loop,  and  the  direction  of  motion 
of  the  ion  are  recorded.  This  is  repeated  until  sufficient  statistical 
significance  is  achieved.  Because  the  ions  do  not  escape  (i.e., 
they  either  interact  or  thermalize),  the  calculations  are  “thick 
target.”  For  deexcitation  line  reactions,  the  resulting  interacting- 
particle  angular  distribution  can  then  be  used  with  the  7-ray  line 
production  code  to  calculate  line  shapes  and  shifts.  The  attenu¬ 
ation  of  the  7-ray  along  the  line  of  sight  is  calculated  from  the 
location  of  the  interaction  along  the  loop.  For  neutron-producing 
reactions,  the  Monte  Carlo  simulation  follows  the  neutron  until  it 
escapes  into  space,  is  captured  on  H  to  contribute  the  2.223  MeV 
capture  line,  reacts  with  3He  without  production  of  radiation  (see 
§  2.4.3),  or  decays.  If  a  neutron-capture  7-ray  is  produced,  its  at¬ 
tenuation  along  the  line  of  site  is  calculated. 

2.2.  Summary  of  the  Acceleration  and  Physical  Parameters 

Table  1  lists  the  acceleration  and  physical  parameters  of  the 
loop  transport  and  interaction  model  that  affect  the  measurable 
quantities  addressed  by  our  calculations.  We  discuss  and  review 
the  range  of  parameter  values  obtained  by  past  analyses. 

2.2.1.  Acceleration  Parameters 

A  fundamental  property  of  accelerated  ions  is  their  energy 
spectrum.  For  the  calculations  performed  here,  we  assume  that 
the  spectral  form  of  the  accelerated  ions  is  an  unbroken  power 
law,  with  all  ion  species  having  the  same  index  (here  denoted  by 
s),  normalized  to  one  proton  with  energy  greater  than  30  MeV. 
Gamma-ray  and  neutron  observations  provide  various  measures 
of  the  power-law  index.  Ramaty  et  al.  (1996)  derived  power-law 
indexes  for  several  flares  observed  with  the  Gamma  Ray  Spec¬ 
trometer  on  the  Solar  Maximum  Mission  ( SMM )  using  two 
such  measures,  the  fluence  ratio  of  the  6.129  MeV  160  to  the 
1.634  MeV  20Ne  lines  and  the  fluence  ratio  of  the  2.223  MeV 
neutron  capture  to  the  4.438  MeV  12 C  lines.  For  each  flare,  he 
found  that  the  indexes  derived  with  the  two  tests  were  the  same, 
within  their  uncertainties  (see  the  Ne/O  discussion  in  §  2.2.2).  We 
will  see  in  §  2.4.4  that  the  first  measure  is  relevant  to  accelerated 
ions  with  energies  from  a  few  to  about  10  MeV  nucleon-1,  while 
the  second  is  relevant  to  ions  with  energies  from  a  few  to  as  much 
as  1 00  MeV  nucleon- 1 ,  implying  that  the  spectmm  is  an  unbroken 
power  law  through  this  entire  energy  range. 
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TABLE  2 

Ambient  and  Accelerated-Ion  Compositions 


Element 

Ambient 

Accelerated 

H . 

1.0 

1.0 

T  lc . 

3.0  x  10-5 

0.5 

;I  le . 

0.1 

0.5 

C . 

2.96  x  10“4 

4.65  x  10  "3 

N . 

7.90  x  10“5 

1.24  x  10  "3 

O . 

6.37  x  10“4 

1.00  x  10  "2 

Ne . 

1.59  x  10“4 

4.55  x  10  "3 

Mg . 

1.25  x  10“4 

5.89  x  10-3 

Al . 

1.00  x  10“5 

1.57  x  10-4 

Si . 

9.68  x  10“5 

4.55  x  10  "3 

S . 

2.03  x  10“5 

9.56  x  10-4 

Ca . 

6.75  x  10“6 

1.06  x  10“4 

Fe . 

8.54  x  10“5 

1.34  x  10  "2 

In  addition  to  using  these  spectral  measures  in  their  analysis  of 
OSSE  data  from  the  1991  June  11  solar  flare,  Murphy  &  Share 
(2000)  used  an  additional  measure:  the  fluence  ratio  of  the 
0.5 1 1  MeV  positron  annihilation  to  the  4.438  MeV  12  C  lines.  We 
will  show  in  §  2.4.4  that  this  measure  is  relevant  to  accelerated 
ions  energies  from  a  few  to  a  few  tens  of  MeV  nucleon-1;  i.e., 
between  the  energy  ranges  of  the  other  two  tests.  They  showed 
that  there  were  time  periods  during  the  June  11  flare  when  the 
spectrum  appeared  to  be  an  unbroken  power  law  throughout  the 
entire  relevant  energy  range,  but  there  were  also  periods  when 
it  was  not,  exhibiting  a  soft  low-energy  component  and  a  hard 
high-energy  component.  Using  similar  spectral  index  measures 
as  those  above,  we  will  show  in  this  paper  that  the  accelerated  ion 
spectrum  for  the  June  4  flare  is  again  consistent  with  an  unbroken 
power  law  up  to  ^100  MeV  nucleon-1  but,  with  the  addition  of 
neutron  measurements  also  available  for  this  flare  (which  are 
relevant  to  ion  energies  up  to  several  hundred  MeV  nucleon- 1 ;  see 
§  2.4.4),  the  spectrum  must  steepen  above  100  MeV  nucleon-1. 

In  our  calculations,  we  will  vary  the  power-law  spectral  index 
from  2  to  6.  In  their  derivation  of  spectral  indexes  from  19  flares 
observed  by  SMM,  Ramaty  et  al.  (1996)  found  most  indexes  fell 
within  this  range,  with  an  average  of  ~4.2.  Spectral  indexes  de¬ 
rived  using  recent  data  from  three  solar  flares  obtained  with  the 
RHESSI  instrument  have  been  harder,  with  indexes  nearer  to  3  or 
less. 

The  accelerated-ion  composition  used  in  the  calculations  is 
summarized  in  Table  2.  We  have  assumed  the  “impulsive  flare” 
accelerated-ion  abundances  defined  by  Ramaty  et  al.  (1996), 
which  reflects  average  heavy-element  abundance  enhancements 
found  in  impulsive  SEP  events  in  space  (Reames  et  al.  1994). 
Accelerated  ion  compositions  derived  from  7-ray  line  mea¬ 
surements  (Murphy  et  al.  1991;  Share  &  Murphy  1998,  1999; 
Mandzhavidze  et  al.  1999)  have  suggested  that  they  are  similar 
to  those  of  such  impulsive  SEP  events.  The  impulsive  compo¬ 
sition  is  coronal  (Reames  1995)  forC,  N,  Ne,  Mg,  Al,  Si,  S,  Ca, 
and  Fe  relative  to  O  but  with  Ne/O,  Mg/O,  Si/O,  and  S/O  in¬ 
creased  by  a  factor  of  3  and  Fe/O  increased  by  a  factor  of  10. 
We  will  assume  two  values  for  a/proton:  0.1  and  0.5.  For  each, 
we  will  maintain  a/O  =  50  and  3He/4He  =  1. 

2.2.2.  Physical  Parameters 

The  ambient  composition  used  in  the  calculations  is  also 
summarized  in  Table  2.  We  have  assumed  coronal  abundances 
(Reames  1995)  for  C,  N,  O,  Mg,  Al,  Si,  S,  Ca,  and  Fe  relative  to 
H.  Abundance  determinations  from  analyses  of  7-ray  lines  for 
several  flares  (Share  &  Murphy  1995;  Ramaty  et  al.  1996)  have 


shown  that  the  composition  of  the  7-ray  line  production  region  is 
similar  to  coronal.  These  are  the  same  ambient  abundances  used 
by  Ramaty  et  al.  (1996),  except  we  assume  4He/H  =  0.10  rather 
than  the  coronal  value  of 0.037.  Mandzhavidze  et  al.  (1999)  have 
shown  that  in  some  flares  the  ambient  4He/H  abundance  ratio 
could  be  >0. 1 .  Grevesse  et  al.  ( 1 996)  list  the  photospheric  4He/ H 
as  0.097. 

We  have  also  taken  Ne/O  =  0.25.  In  their  analysis  of  19  SMM 
7-ray  line  flares,  Ramaty  et  al.  ( 1 996)  showed  that  to  maintain  an 
unbroken  power-law  spectrum  from  ~2  to  ~50  MeV  nucleon-1, 
the  ambient  Ne/O  abundance  ratio  had  to  be  ~0.25,  higher  than 
the  accepted  photospheric  value  of  0.15.  Similar  Ne/O  enhance¬ 
ments  have  been  seen  in  other  measurements  of  flare  data;  for 
example,  Schmelz  (1993)  found  Ne/O  of  0.32  using  observations 
of  a  solar  flare  obtained  with  the  Flat  Crystal  Spectrometer  on 
SMM.  Drake  &  Testa  (2005)  found  Ne/O  =  0.41  for  a  sample  of 
21  stars  located  within  100  pc  of  the  Sun  that  had  been  observed 
with  the  High  Energy  Transmission  Grating  on  the  Chandra 
X-Ray  Observatory.  But  the  Ne/O  abundance  ratio  derived  in  a 
recent  analysis  of  FCS  data  that  excluded  data  taken  during  flares 
and  long-duration  events  (Schmelz  et  al.  2005)  was  consistent 
with  the  accepted  0.15  value.  Schmelz  et  al.  (2005)  also  reconsid¬ 
ered  older  data  obtained  from  full-Sun  observations  and  again 
found  Ne/O  consistent  with  the  accepted  value.  Considering  that 
the  Ne/O  enhancements  could  be  flare-related,  Shemi  (1991) 
suggested  that  preflare  soft  X-ray  radiation  preferentially  ionizes 
Ne,  which  could  then  be  transported  to  higher  levels  of  the  atmo¬ 
sphere.  Since  the  Chandra  spectra  analyzed  by  Drake  &  Testa 
(2005)  are  most  likely  dominated  by  photons  from  very  intense 
flares,  their  enhanced  neon  abundances  might  result  from  a  sim¬ 
ilar  mechanism.  However,  Schmelz  et  al.  (2005)  point  out  that 
this  is  not  likely  since  low-FlP  elements  in  active  stellar  coronae 
are  generally  depleted  rather  than  enhanced,  as  they  are  in  the  so¬ 
lar  corona.  Here  we  will  use  the  0.25  value  derived  from  the  solar- 
flare  7-ray  measurements.  In  our  analysis  of  the  June  4  flare  7-ray 
line  data  (see  §  3. 1),  we  will  see  that  the  assumption  of  an  unbroken 
power  law  remains  valid  through  the  accelerated  ion  energy  range 
relevant  to  those  measurements  (^2-50  MeV  nucleon-1).  The 
addition  of  the  neutron  data  will  show  that  the  spectrum  could 
steepen  above  those  energies,  however  (see  §  3.1). 

We  will  investigate  the  effect  on  the  measurable  quantities  for 
two  models  of  the  solar  atmosphere:  the  sunspot  region  atmo¬ 
sphere  of  Avrett  (1981)  and  a  hybrid  model  that  could  be  more 
similar  to  the  atmosphere  during  an  active  period.  In  the  sunspot 
model  of  Avrett,  the  density  at  heights  greater  than  several  hun¬ 
dred  kilometers  is  more  tenuous  than  that  of  the  quiet-Sun  model 
of  Vemazza  et  al.  (1981)  but  below  this  height  it  is  greater.  Dur¬ 
ing  a  flare,  the  atmosphere  could  be  considerably  different  still 
(Share  et  al.  2004).  Aschwanden  et  al.  (2002)  have  used  RHESSI 
measurements  of  electron  bremsstrahhmg  to  estimate  the  density- 
height  profile  in  a  solar  flare  and  found  that  the  density  is  enhanced 
relative  to  quiet-Sun  models  at  heights  greater  than  ~  1 000  km  and 
the  density  scale  height  is  considerably  larger.  We  construct  a  hy¬ 
brid  solar-flare  model  (the  “RHESSI/ Avrett  model”)  consisting  of 
the  Avrett  density  profile  below  1000  km  and  the  Aschwanden 
profile  above.  We  use  the  Avrett  sunspot  model  at  photospheric 
heights  since  flare  loop  structures  are  generally  found  above  sun¬ 
spots  (e.g.,  Hurford  et  al.  2003;  Schrijver  et  al.  2006).  The  en¬ 
hanced  chromospheric  density  during  flares  suggested  by  the 
RHESSI  measurements  could  be  due  to  filling  of  the  loop  by  ac¬ 
tive  processes  prior  to  the  main  impulsive  event.  The  two  mod¬ 
els,  along  with  the  Vemazza  model,  are  shown  in  Figure  1 ,  where 
the  larger  chromospheric  density  and  scale  height  of  the  RHESSI/ 
Avrett  model  is  clearly  seen.  This  is  the  first  time  that  the 
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Fig.  1 . — Density  vs.  height  in  the  solar  atmosphere  for  the  three  atmospheric 
models  considered  here:  Vemazza  (  Vemazza  et  al.  1981),  Avrett  (Avrett  1981), 
and  RHESSI/ Avrett. 


dependence  of  solar-flare  7-ray  and  neutron  production  on  the 
atmospheric  model  has  been  systematically  explored,  thus  al¬ 
lowing  the  possibility  to  learn  about  the  structure  of  the  solar 
atmosphere  during  flares  from  high-energy  measurements. 

The  parameters  directly  associated  with  ion  transport  in  the 
loop  model  are  the  level  of  PAS,  the  degree  of  loop  magnetic 
convergence,  and  loop  length.  We  will  vary  the  level  of  PAS  from 
none  (A  — >  00)  to  saturated.  At  saturation  (A  ss  20),  further  de¬ 
crease  of  A  has  no  effect  on  the  temporal  or  angular  dependence  of 
the  interactions  (see  discussion  below).  We  will  vary  the  magnetic 
convergence  from  none  ( 8  =  0)  to  very  strong  (6  =  0.45).  For  the 
Avrett  (1981)  atmosphere,  b  =  0.1  and  0.2  correspond  to  an  in¬ 
crease  of  the  magnetic  field  from  the  corona  to  the  photosphere 
by  factors  of  ~5  and  ~20,  respectively.  The  loop  length  will  be 
varied  by  an  order  of  magnitude  from  1 1 ,500  to  1 1 5,000  km.  Es¬ 
timates  of  loop  lengths  based  on  hard  X-ray  timing  and  footpoint 
imaging  (e.g.,  Aschwanden  et  al.  1996;  White  et  al.  2002)  range 
from  ~  10,000  to  80,000  km. 

Nuclear  reactions  and  the  yields  of  the  various  products  in 
solar  flares  are  not  influenced  by  the  temperature  of  the  inter¬ 
action  region.  At  the  temperatures  expected  in  solar  flares,  any 
deexcitation-line  Doppler  broadening  is  dominated  by  the  recoil 
velocities  of  the  excited  nuclei  resulting  from  the  impact  of  the 
nonthermal  accelerated  particles.  For  the  4.438  MeV  12C  de¬ 
excitation  line,  the  FWHM  due  to  recoil  is  ~80  keV,  while  even  at 
1 06  K  the  thermal  FWHM  is  <  1  keV.  The  width  of  the  2.223  keV 
neutron-capture  line  does  reflect  the  temperature  of  the  region 
where  the  neutrons  are  captured  since  the  neutrons  are  typically 
thermalized  before  capture,  but  since  the  captures  occur  deep  in 
the  photosphere  where  the  temperature  is  <  1 04  K,  the  FWHM  is 
~10  eV,  much  less  than  typical  detector  energy  resolution  at 
these  energies.  The  yield  of  secondary  positrons  from  production 
of  radioactive  nuclei  and  pions  is  similarly  not  influenced  by  the 
temperature,  but  the  yield  of  the  subsequent  annihilation  pho¬ 
tons  and  the  shape  of  the  annihilation  line  is.  Positron  annihila¬ 
tion  was  thoroughly  studied  by  Murphy  et  al.  (2005). 

2.2.3.  Nuclear  Interaction  Cross  Sections 

The  development  of  the  theory  of  7-ray  line  production  mo¬ 
tivated  a  systematic  effort  of  laboratory  measurements  of  pro¬ 
duction  cross  sections  for  the  most  important  lines  in  solar  flares 
with  resulting  estimated  uncertainties  of  less  than  20%  (e.g., 
Dyeretal.  1981, 1985).  The  7-ray  production  code  developed  by 


Ramaty  et  al.  (1979)  now  includes  over  100  individual  lines.  The 
cross  sections  for  all  of  these  lines  have  recently  been  re-evaluated 
and  updated  with  new  laboratory  measurements  (  Kozlovsky  et  al. 
2002).  The  cross  sections  are  generally  strongly  peaked  in  the  en¬ 
ergy  range  10-30  MeV  nucleon-1. 

The  7-ray  line  production  code  includes  accurate  kinematics 
to  calculate  the  direction  of  the  recoil  excited  nucleus  in  the  center- 
of-mass  frame  and  the  angular  distribution  ofthe  7-ray  emission  in 
the  frame  of  the  excited  nucleus  (see  Ramaty  et  al.  1979  and  the 
more  recent  work  of  Kiener  et  al.  2001).  Anisotropic  emission  in 
the  excited  nucleus  rest  frame  generally  will  not  affect  line  fluence 
detemiinations  unless  the  accelerated  particles  are  truly  beamed. 
In  addition,  the  velocities  of  the  excited  nuclei  are  generally  too 
small  to  result  in  significant  anisotropy  of  the  emission  in  the 
frame  of  the  Sun  due  to  relativistic  beaming.  However,  the  veloc¬ 
ities  are  large  enough  to  produce  measurable  line  Doppler  shifts 
when  the  recoil  angular  distribution  is  sufficiently  anisotropic. 

In  addition  to  the  explicit  7-ray  lines  included  in  our  deexcita¬ 
tion  line  code,  nuclear  reactions  involving  heavier  nuclei  also 
produce  a  large  number  of  weaker,  closely  spaced  lines  and  a 
true  nuclear  continuum.  This  additional  emission  is  currently  in¬ 
cluded  in  the  line  production  code  according  to  the  prescription 
given  by  Ramaty  et  al.  (1979).  We  are  currently  improving  this 
modeling  using  both  measurements  and  calculations  of  this  emis¬ 
sion.  Accurate  modeling  of  this  emission  is  not  critical  for  the  pur¬ 
poses  of  this  paper  since  we  will  focus  only  on  measurements  of 
the  strongest  narrow  lines. 

Hua  et  al.  (2002)  thoroughly  reevaluated  and  revised  the  neu¬ 
tron  production  cross  sections  and  kinematics.  With  these  new 
approximations  for  the  angular-  and  energy-dependent  cross  sec¬ 
tions,  the  calculated  neutron  yields  are  in  much  better  agreement 
with  laboratory  measurements.  The  estimated  uncertainty  in  the 
calculated  neutron  yield  over  the  neutron  energy  range  relevant 
to  solar  flare  measurements  is  ~30%. 

2.3.  Effect  of  the  Acceleration  and  Physical  Parameters 
on  the  Nuclear  Interactions 

Except  for  neutrons  that  escape  from  the  Sun,  the  products  of 
nuclear  reactions  (e.g.,  excited  and  radioactive  nuclei,  pions,  and 
neutrons)  are  not  directly  observable  but  reveal  themselves  through 
secondary  processes  producing  emissions  that  are  observable  (e.g., 
deexcitation  7-rays,  annihilation  7-rays  after  positron  emission, 
pion-decay  radiation,  and  neutron  capture  7-rays).  The  proper¬ 
ties  of  an  interacting  accelerated  ion  that  most  significantly  affect 
the  measurable  quantities  associated  with  the  observable  emis¬ 
sions  are  (1)  the  time  of  the  interaction  after  release  from  the 
acceleration  region,  (2)  the  ion  direction  of  motion  at  interaction, 
(3)  where  in  the  loop  the  interaction  occurs,  and  (4)  the  ion 
kinetic  energy  at  interaction.  The  parameters  listed  in  Table  1 
combine  to  affect  all  of  these  properties  to  varying  degrees.  In  the 
following  subsections  we  discuss  how  these  interacting  ion  prop¬ 
erties  depend  on  the  parameters,  using  calculations  of  the  pro¬ 
duction  of  the  4.438  MeV  excited  state  of  12C  for  illustration.  The 
results  are  qualitatively  valid  for  other  deexcitation  lines  and  for 
neutron  production.  The  total  number  of  possible  parameter  com¬ 
binations  is  too  large  to  fully  explore  so  the  parameter  ranges  and 
values  discussed  are  chosen  to  effectively  convey  the  most  impor¬ 
tant  dependences.  We  discuss  the  interaction  time  history  in  §  2.3 . 1 , 
the  interacting  ion  angular  distribution  in  §  2.3.2,  and  the  distri¬ 
bution  of  the  interactions  in  the  loop  in  §  2.3.3.  We  defer  the 
discussion  of  the  energies  of  the  interacting  ions  to  §  2.4.4  after 
discussing  the  effect  of  the  parameters  on  the  various  observable 
emissions  resulting  from  the  nuclear  interactions  in  §§  2.4. 1 , 2.4.2, 
and  2.4.3. 
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Fig.  2. — Dependence  of  the  nuclear  interaction  rate  time  history  on  various  parameters:  (a)  magnetic  convergence  (6);  (b)  PAS  (2);  (c)  loop  length  (. L );  (d)  spectral 
index  (5);  and  ( e )  atmospheric  model.  Panel  (/)  shows  the  range  of  the  effect  that  PAS  has  on  the  Avrett  and  RHESSI/ Avrett  atmospheric  models.  For  each  panel,  the  other 
parameters  are  fixed  at  the  values  noted  in  the  text.  All  yields  are  normalized  to  1  accelerated  proton  of  energy  greater  than  30  MeV. 


In  the  following  sections,  figures  are  presented  showing  the 
dependence  of  the  interacting  ion  properties  on  a  particular  ac¬ 
celeration  or  physical  parameter.  Except  for  the  parameter  being 
varied,  and  unless  otherwise  indicated,  the  other  parameters  will 
be  held  fixed  at  a  set  of  standard  values:  a/proton  =  0.5,  s  =  4, 
X  =  300,  8  =  0.2,  L  =  11 ,500  km,  and  the  atmosphere  of  Avrett 
(1981). 

2.3.1.  Nuclear  Interaction  Rate  Time  History 

The  time  history  of  the  nuclear  interaction  rate  is  governed  by 
both  the  time  history  of  the  ion  acceleration  release  rate  (aion(0) 
and  transport  in  the  loop.  Here  we  will  only  consider  the  effects 
of  transport,  assuming  instantaneous  release  of  all  particles  at 
t  =  0.  (When  the  June  4  flare  data  are  analyzed  in  §  3,  we  will 
discuss  estimates  for  aion(t).)  In  the  absence  of  magnetic  conver¬ 
gence  (5  =  0),  ions  do  not  mirror  and  the  nuclear  interaction  time 
history  depends  only  on  the  energy  loss  rate  in  the  lower  chro¬ 
mosphere  and  upper  photosphere  where  the  density  is  greatest 
and  most  of  the  interactions  occur.  In  the  presence  of  magnetic 
convergence  ( 8  ^  0),  ions  with  initial  pitch  angles  lying  close  to 
the  loop  axis  also  do  not  mirror  (they  are  within  the  “loss  cone”) 


and  behave  similarly.  Ions  outside  the  loss  cone  mirror  and  lose 
energy  much  more  slowly  as  they  traverse  the  low-density  co¬ 
rona.  These  interactions  thus  occur  on  longer  timescales,  with 
the  timescale  increasing  with  increasing  convergence.  Figure  2 a 
shows  12 C  line  production  time  histories  for  8  =  0.0,  0.20,  and 
0.45.  For  the  time  histories  with  8  ^  0,  the  excesses  at  times  less 
than  1  s  are  due  to  ions  within  the  loss  cone  interacting  on  their 
first  pass  down  the  loop  leg.  For  all  of  these  calculations,  we  have 
assumed  no  PAS;  i.e.,  X  — >  oo. 

Pitch  angle  scattering  causes  the  loss  cone  to  be  continuously 
repopulated.  As  a  result,  the  nuclear  interaction  rate  increases 
at  early  times  and  correspondingly  decreases  at  later  times.  Fig¬ 
ure  2b  shows  12C  line  production  time  histories  for  PAS  levels  of 
none  {X  — >  oo),  intermediate  (A  =  300),  and  high  (X  =  20).  The 
time  history  is  no  longer  affected  by  decreasing  X  when  the  rate 
of  loss-cone  replenishment  exceeds  the  rate  of  nuclear  reactions 
in  the  loss  cone.  This  saturation  effect  occurs  for  X  less  than  ~20. 

Increasing  the  loop  length  increases  the  timescale  of  the  inter¬ 
action  rate  since  mirroring  particles  spend  more  time  at  lower 
coronal  densities  where  nuclear  reactions  are  less  likely.  Figure  2c 
shows  the  effect  on  the  interaction  time  history  due  to  loop  length. 
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The  delays  of  the  excesses  at  ~  1  s  and  less  ( for  ions  within  the  loss 
cone  interacting  on  their  first  pass  down  the  loop  leg)  are  directly 
due  to  the  longer  transit  times  resulting  from  the  longer  loop 
lengths. 

The  accelerated  ion  spectral  index  has  no  significant  effect  on 
the  interaction  rate  time  history  at  times  later  than  about  1  s  after 
particle  release.  However,  in  harder  spectra  there  is  a  larger  frac¬ 
tion  of  ions  with  higher  energies  producing  the  reactions  (see 
§  2.4.4).  Before  about  1  s,  this  results  in  ions  arriving  earlier  at 
the  footpoints  on  their  first  transport  down  the  loop  leg.  This  can 
be  seen  in  Figure  2d,  where  time  histories  for  5  =  3,4,  and  5  are 
shown.  Similarly,  different  energy  dependences  of  the  interac¬ 
tion  cross  sections  affect  the  time  history  only  before  about  1  s. 

The  differences  of  the  Avrett  and  Vemazza  atmospheric  mod¬ 
els  have  negligible  effect  on  the  interaction  time  history.  Only 
when  the  model  is  significantly  different,  such  as  the  RHESSII 
Avrett  model,  is  the  time  history  affected,  with  a  shortened  time- 
scale  due  to  the  dramatic  chromospheric  and  lower  corona  den¬ 
sity  enhancement  as  shown  in  Figure  2e.  Another  aspect  of  the 
time  history  dependence  on  the  atmospheric  model  is  shown  in 
Figure  If.  At  times  less  than  a  few  seconds,  the  extent  of  the 
variation  of  the  time  history  as  PAS  is  varied  from  none  to  sat¬ 
urated  is  significantly  less  for  the  RHESSU Avrett  model  than  it  is 
for  the  Avrett  model.  A  similar  reduced  extent  of  the  variation  is 
found  as  8  is  varied.  This  will  impact  the  ability  to  distinguish 
among  possible  values  of  8  and  X  when  data  from  the  June  4  flare 
are  analyzed  in  §  3 . 1  under  the  assumption  of  the  RHESSU Avrett 
model. 

2.3.2.  Interacting  Ion  Angular  Distribution 

As  mentioned  in  §  2.1,  calculations  based  on  the  loop  model 
produce  angular  distributions  for  the  interacting  ions  that  can  be 
related  to  physical  parameters,  in  contrast  to  the  ad  hoc  distribu¬ 
tions  assumed  previously.  In  the  absence  of  magnetic  conver¬ 
gence  (8  =  0),  there  is  no  mirroring  and  the  angular  distribution 
of  the  interacting  ions  (initially  isotropic  when  released  at  the  top 
of  the  loop)  will  be  downward  isotropic  regardless  of  level  of  PAS 
or  loop  size.  In  the  presence  of  magnetic  convergence  (8  f  0)  but 
no  PAS  (X  —>  oo),  ions  with  pitch  angles  too  small  to  mirror  enter 
the  “loss  cone”  and  either  undergo  a  nuclear  reaction  as  they  are 
moving  downward  or  are  thermalized.  But  most  ions  have  larger 
initial  pitch  angles,  which  results  in  mirroring,  and  the  interactions 
occur  mostly  at  the  ion  mirror  points  where  the  density  is  greatest. 
Thus,  most  ions  are  moving  parallel  to  the  solar  surface  when  they 
interact  (i.e.,  a  “fan  beam”  angular  distribution;  see,  e.g.,  Hua  & 
Lingenfelter  1987b)  with  a  relatively  small  contribution  from 
downward-moving  ions.  The  downward  contribution  becomes 
even  less  important  as  8  increases.  Figure  3 a  shows  the  effect  of 
magnetic  convergence  on  the  interacting  ion  angular  distribu¬ 
tion  for  8  =  0.0,  0.20,  and  0.45.  For  all  of  these  calculations, 
there  is  no  PAS  {X  — *  oo). 

By  continuously  repopulating  the  loss  cone,  the  effect  of  PAS 
is  to  produce  more  interactions  involving  downward-directed 
ions.  This  is  seen  in  Figure  3 b,  where  angular  distributions  are 
shown  for  the  same  convergences  as  Figure  2 a  ( 8  =  0.0,  0.20, 
and  0.45)  but  for  X  =  300.  As  PAS  is  increased  ( X  decreased), 
more  ions  interact  in  the  loss  cone  and  the  distribution  becomes 
more  downward-directed  until  saturation  is  reached  ( X  ~  20). 
This  is  shown  in  Figure  3c  for  X  — >  oo,  X  =  300,  and  X  =  20. 

The  angular  distribution  of  the  accelerated  ions  when  they 
interact  is  not  significantly  affected  by  the  spectral  index  or  the 
energy  dependence  of  the  interaction  cross  section.  The  distribu¬ 
tions  produced  by  the  three  atmospheric  models  considered  here 
(Vemazza,  Avrett,  and  RHESSU Avrett)  were  also  very  similar, 


Fig.  3. — Dependence  of  the  interacting  particle  angular  distribution  on  vari¬ 
ous  parameters:  (a)  magnetic  convergence  (6)  with  no  PAS  (2  — >  oo);  ( b )  mag¬ 
netic  convergence  (6)  with  X  =  300;  and  (c)  PAS  (with  6  =  0.20).  For  each  panel, 
the  other  parameters  are  fixed  at  the  values  noted  in  the  text,  [l  —  cos  (9),  where  9  is 
the  angle  between  the  normal  to  the  solar  surface  and  the  direction  of  the  ion.  9  —  0° 
(p  —  1)  is  directed  outward  from  the  Sun,  9  =  180°  (/i  =  —  1)  is  directed  inward 
into  the  Sun,  and  9  =  90°  (/i  =  0)  is  directed  parallel  to  the  surface  of  the  Sun.  The 
distributions  are  normalized  to  unity. 

except  for  low  levels  of  PAS  (large  X),  where  the  RHESSU Avrett 
distribution  had  a  slightly  enhanced  downward-directed  compo¬ 
nent  resulting  from  enhanced  interactions  of  particles  within  the 
loss  cone. 

2.3.3.  Interaction  Depth  Distribution 

In  the  absence  of  magnetic  convergence  ( 8  =  0),  there  is  no 
mirroring  and  the  depth  distribution  directly  reflects  the  column 
density  required  for  the  accelerated  particles  to  interact  as  they 
move  downward  through  the  solar  atmosphere.  In  the  presence 
of  magnetic  convergence  ( 8  f  0),  mirroring  results  in  interac¬ 
tions  occurring  at  higher  elevations  (and  therefore  lower  densities) 
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Fig.  4. — Dependence  of  the  depth  distribution  of  the  nuclear  interactions  on  various  parameters:  ( a )  magnetic  convergence  (<5);  ( b )  PAS  (A);  (c)  spectral  index  (s); 
(i d )  atmospheric  model;  ( e )  interaction  cross  section;  and  (/)  accelerated  cn/proton  ratio.  For  each  panel,  the  abscissa  is  the  density  in  the  solar  atmosphere,  which  can 
be  related  to  depth  through  the  atmospheric  model.  Also  for  each  panel,  the  other  parameters  are  fixed  at  the  values  noted  in  the  text. 


as  more  and  more  particles  are  prevented  from  penetrating  the 
lower  atmosphere.  This  is  seen  in  Figure  4 a,  where  the  fraction  of 
interactions  occurring  deeper  than  a  given  density  is  plotted  versus 
density  in  the  solar  atmosphere  for  6  =  0.0,  0.20,  and  0.45. 

Even  with  minimal  convergence,  the  lack  of  PAS  (2  — >  oo)  re¬ 
sults  in  a  significant  fraction  of  the  interactions  occurring  at  low 
densities.  This  can  be  seen  in  Figure  4b,  where  depth  distribu¬ 
tions  for  no  PAS  (2  — >  oo),  2  =  300,  and  saturated  PAS  (2  =  20) 
are  shown.  The  curve  for  no  PAS  shows  that  almost  20%  of 
the  interactions  occur  at  densities  less  than  10 10  H  cnr3.  As 
PAS  is  increased  (2  decreased),  the  bulk  of  the  interactions  move 
to  higher  densities  as  more  interactions  occur  within  the  loss 
cone. 

Because  the  nuclear  interactions  occur  predominantly  at  chro¬ 
mospheric  depths  or  greater,  transport  through  the  low-density 
corona  has  little  effect  on  where  the  interactions  will  ultimately 
occur.  As  a  result,  the  effect  of  loop  length  on  the  interaction  depth 
distribution  is  minimal,  as  long  as  the  other  parameters  are  not 
changed.  Specifically,  to  maintain  the  relative  level  of  PAS  in  the 
calculations  when  the  loop  length  is  increased,  the  physical  scat¬ 
tering  mean  free  path  (A)  must  be  increased  accordingly.  For  a 


fixed  A,  however,  increasing  the  loop  length  results  in  decreasing 
2,  the  effect  of  which  was  discussed  above. 

Because  higher  energy  ions  tend  to  interact  farther  along  their 
paths,  harder  particle  spectra  result  in  interactions  occurring  at 
higher  densities  since  more  ions  of  higher  energies  are  producing 
the  interactions  (see  §  2.4.4).  This  is  seen  in  Figure  4c  for  s  =  3, 
4,  and  5. 

When  the  depth  distribution  is  considered  as  a  function  of 
atmospheric  density,  the  effect  of  different  atmospheric  density 
profiles  is  minimal  as  seen  in  Figure  4 d,  where  the  depth  distri¬ 
bution  is  shown  for  the  three  models  considered  here  ( Vemazza, 
Avrett,  and  RHESSI/AwreXX).  This  is  due  partly  to  the  fact  that  for 
convenience  we  have  parameterized  the  magnetic  convergence 
in  terms  of  the  pressure  so  that,  for  a  given  6,  the  mirror  point 
tends  to  remain  near  the  same  density  as  the  atmospheric  model 
is  changed.  In  terms  of  physical  depth,  however,  the  distribution 
does  change  significantly  as  the  model  is  changed,  as  seen  in  Fig¬ 
ure  5,  where  the  distribution  as  a  function  of  depth  in  the  atmo¬ 
sphere  is  shown  for  the  three  model  atmospheres. 

Interactions  whose  cross  sections  have  high  thresholds  tend  to 
interact  at  higher  densities,  again  the  consequence  of  the  longer 
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Fig.  5. — Dependence  of  the  depth  distribution  of  the  nuclear  interactions  on 
the  atmospheric  model  as  in  Fig.  4 d  but  with  the  abscissa  being  depth  in  the 
atmosphere.  For  each  panel,  the  other  parameters  are  fixed  at  the  values  noted 
in  the  text. 


ranges  of  higher  energy  ions.  This  is  seen  in  Figure  4e,  where  the 
depth  distributions  for  interactions  with  two  different  cross  section 
energies  are  shown:  a  low  threshold  energy  (~1  MeV  nucleon-1, 
similar  to  that  of  the  1 .634  MeV  20Ne  deexcitation  line)  and  a  high 
threshold  energy  (~25  MeV  nucleon-1,  similar  to  those  of  neutron- 
producing  p-a  reactions). 

This  dependence  of  the  depth  distribution  on  the  interaction 
cross  section  also  explains  why  the  distribution  is  affected  by  the 
accelerated  a/proton  ratio,  as  shown  in  Figure  4 f  This  depen¬ 
dence  results  from  the  different  threshold  energies  of  the  a  and 
proton  cross  sections  for  producing  the  4.438  MeV  excited  state 
of  12C.  The  a-particle  interactions  have  a  lower  threshold  and,  as 
just  discussed,  such  interactions  occur  at  lower  densities  (see 
Fig.  4e).  When  the  a/proton  ratio  is  high,  a  larger  fraction  of  the 
line  yield  is  due  to  such  reactions  and  the  depth  distribution  shifts 
to  lower  densities. 

2.4.  Effect  of  the  Acceleration  and  Physical  Parameters 
on  the  Observable  Quantities 

In  §  2.3  we  discussed  how  the  various  acceleration  and  phys¬ 
ical  parameters  affect  those  aspects  of  the  nuclear  interactions 
that  most  significantly  affect  the  subsequent  observable  emissions. 
In  the  next  three  subsections  we  discuss  how  the  parameters  affect 
the  measurable  quantities  associated  with  deexcitation  lines,  es¬ 
caping  neutrons,  and  the  neutron  capture  line,  respectively.  In  a 
final  subsection  we  discuss  what  accelerated  ion  energies  are  re¬ 
sponsible  for  each  of  these  emissions. 

2.4.1.  Deexcitation  Line  Production 

In  this  section  we  discuss  the  measurable  quantities  associated 
with  “narrow”  nuclear  deexcitation  lines  produced  by  accelerated 
proton  and  a-particle  interactions  with  ambient  C  and  heavier  nu¬ 
clei.  For  such  deexcitation  lines,  the  measurable  quantities  are  the 
line  fluence,  shape,  central  energy,  and  emission  time  history. 
How  the  acceleration  and  physical  parameters  affect  each  of 
these  quantities  is  discussed. 

2. 4. 1.1.  Deexcitation  Line  Yields  and  Yield  Ratios 

The  only  physical  parameter  significantly  affecting  time- 
integrated  narrow  deexcitation  line  yields  (and  also,  therefore, 
line  ratios)  is  the  relative  abundances  of  the  ambient  elements 
heavier  than  He.  Yields  are  not  affected  by  the  physical  parameters 


Fig.  6. — Calculated  yields  of  the  1.634  MeV  (20Ne),  4.438  MeV  (12C), 
6. 129  MeV  (160),  and  ~0.45  MeV  a-  a  complex  as  functions  of  the  accelerated- 
ion  power-law  spectral  index  for  accelerated  a/proton  —  0.5  (thick  lines)  and 
0. 1  ( thin  lines).  Deexcitation  line  yields  do  not  depend  on  loop  parameters.  All 
yields  are  normalized  to  1  accelerated  proton  of  energy  greater  than  30  MeV. 

associated  with  the  loop  itself,  such  as  length,  PAS  or  magnetic 
convergence.  This  is  because  in  the  loop  model  accelerated  ions 
do  not  escape  from  the  loop,  so  the  interactions  are  “thick  target” 
overall.  In  addition,  the  lines  are  not  produced  sufficiently  deep  in 
the  solar  atmosphere  for  significant  attenuation  to  occur.  Only 
when  the  flare  is  located  at  or  beyond  the  limb  of  the  Sun  is  there 
any  significant  attenuation  of  any  line.  Also,  the  interacting  angu¬ 
lar  distribution  does  not  significantly  affect  line  ratios  because  the 
recoil  velocities  of  the  excited  nuclei  are  nonrelativistic,  so  there 
is  no  significant  beaming  of  the  emission.  The  anisotropic  emis¬ 
sion  of  photons  in  the  rest  frame  of  the  recoil  nucleus  is  signif¬ 
icant  only  for  very  beamed  distributions.  Therefore,  calculated 
total  production  yields  and  yield  ratios  can  be  used  for  flares 
occurring  anywhere  on  the  solar  disk. 

The  acceleration  parameters  that  affect  narrow  deexcitation  line 
yields  (and  ratios)  are  the  accelerated  a/proton  ratio  and  the  spec¬ 
tral  index.  Both  accelerated  protons  and  a-particles  contribute  to 
essentially  all  reactions  producing  excited  nuclei  via  direct  excita¬ 
tion  or  spallation  and  the  relative  contributions  from  the  two  par¬ 
ticles  vary  according  to  the  a  /proton  ratio  and  the  spectral  index 
(through  the  reaction  cross  sections).  The  acceleration  parameters 
that  affect  the  relative  intensities  of  broad  deexcitation  lines  are  the 
accelerated  ion  composition  and  the  spectral  index.  The  only  phys¬ 
ical  parameter  affecting  broad  deexcitation  line  yields  and  line  ra¬ 
tios  is  the  ambient  He/H  ratio. 

Using  the  updated  cross  sections  of  Kozlovsky  et  al.  (2002)  we 
show  in  Figure  6  the  yields  of  the  1.634  MeV  20Ne,  4.438  MeV 
12C,  and  6.129  MeV  160  lines  as  functions  of  the  accelerated-ion 
power-law  spectral  index  for  a/p  =  0.1  ( thin  curves )  and  0.5 
(thick  curves).  These  yields  include  the  direct  excitation  reac¬ 
tions  and  all  spallation  reactions  producing  the  excited  nucleus. 
While  each  of  these  yields  is  dominated  by  the  line  of  the  nu¬ 
cleus  noted,  we  also  include  weaker  lines  at  close-lying  central 
energies.  Specifically,  the  1 .634  MeV  20Ne  yield  also  includes  the 
lines  at  1.600  from  23Mg,  1.635  MeV  from  14N,  and  1.636  MeV 
from  23Na  in  addition  to  the  20Ne  line.  The  4.438  MeV  12C  yield 
also  includes  the  4.444  MeV  line  from  1  *B  in  addition  to  the  12C 
line.  The  6.129  MeV  160  yield  also  includes  the  6.175  MeV 
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Fig.  7. — Calculated  fluence  ratios  of  the  6. 129  MeV  (160)/1 .634  MeV  (20Ne) 
and  2.223  MeV  (neutron-capture)/4.438  MeV  (12C)  lines  as  functions  of  the 
power-law  spectral  index  for  accelerated  a/proton  =  0.5  {thick  lines)  and  0.1 
(thin  lines).  The  latter  ratio  is  for  a  disk-centered  flare  and  is  shown  for  both  no 
PAS  (dashed  line)  and  saturated  PAS  (solid  line).  For  each  level  of  PAS,  the  ratio 
is  shown  for  magnetic  convergences  of  6  =  0.0  and  0.45,  as  indicated  in  the 
figure. 

line  from  150,  the  6.322  MeV  line  from  15N,  and  the  6.337  and 
6.476  MeV  lines  from  nC  in  addition  to  the  160  line.  Including 
these  additional  lines  is  appropriate  when  comparing  these 
calculated  yields  with  yields  measured  with  medium-resolution 
Nal  detectors  such  as  OSSE  and  SMM.  Higher  resolution  de¬ 
tectors  such  as  RHESSI  may  be  able  to  separate  the  individual 
lines.  The  shapes  of  the  curves  in  Figure  6  result  from  (1)  the 
relationship  between  the  accelerated  ion  spectrum  and  the  en¬ 
ergy  dependence  of  the  cross  section  (of  the  main  line  and  also 
of  the  additional  lines  included  in  the  yield)  and  (2)  the  nor¬ 
malization  of  1  proton  above  30  MeV.  The  minimum  exhibited 
by  the  curves  is  an  artifact  due  to  this  normalization  and  the 
index  at  which  the  minimum  occurs  changes  as  the  normaliza¬ 
tion  energy  is  changed.  The  choice  of  this  normalization  energy 
is  arbitrary  and  is  set  at  30  MeV  because  that  is  near  the  energy 
where  many  of  the  deexcitation  line  production  cross  sections 
peak. 

As  first  suggested  by  Lingenfelter  &  Ramaty  (1967),  the  flu¬ 
ence  ratio  of  two  lines  whose  production  cross  sections  have 
markedly  different  energy  dependences  is  sensitive  to  the  steep¬ 
ness  of  the  accelerated-ion  spectrum.  The  6.129  MeV  160  and 
the  1 .634  MeV  20Ne  lines  are  uniquely  useful  in  this  way  because 
(1)  the  threshold  for  the  20Ne  line  is  unusually  low  (~2  MeV  for 
proton  reactions  and  ~1  MeV  nucleon-1  for  a  reactions)  and  is 
very  different  from  that  of  the  160  line  (~7  MeV  for  proton  reac¬ 
tions  and  ~2.5  MeV  nucleon-1  for  a  reactions),  (2)  they  are  both 
strong  lines,  and  (3)  they  both  arise  from  elements  with  high  FIP 
and  are  therefore  expected  to  be  free  of  the  abundance  varia¬ 
tions  associated  with  low  FIP  elements  (Murphy  et  al.  1991; 
Share  &  Murphy  1995;  Ramaty  et  al.  1995).  In  addition  to  the 
time-integrated  measure  of  the  spectral  index,  the  time  depen¬ 
dence  of  the  ratio  also  provides  an  instantaneous  measure  because 
both  of  these  lines  are  prompt  (the  half-lives  of  the  excited  nuclei 


Index 

Fig.  8. — Calculated  fluence  ratio  of  the  cn-a  complex  and  the  sum  of  the 
strongest  high-FIP  lines  (1.63, 4.44,  6.13,  and  ~7  MeV  from20Ne,  12C,  160,  and 
160,  respectively)  as  a  function  of  spectral  index. 

being  <  1 0  10  s).  (Note  that  there  is  no  significant  difference  of  the 
interaction  time  histories  due  to  the  different  energy  dependences 
of  the  cross  sections;  see  §  2.3. 1 .)  Using  this  technique,  power-law 
spectral  indexes  have  been  determined  for  a  number  of  flares 
observed  by  several  detectors  over  the  past  two  solar  cycles  (e.g., 
Ramaty  etal.  1996;  Murphy  etal.  1997;  Ramaty  &  Mandzhavidze 
2000).  This  ratio  was  used  by  Ramaty  &  Mandzhavidze  (2000) 
to  establish  the  near  energy  equivalency  of  accelerated  ions  and 
electrons. 

Figure  7  shows  the  yield  ratio  of  the  6.129  MeV  160  and 
1.634  MeV  20Ne  lines  as  a  function  of  power-law  spectral  index 
for  alp  =  0.1  {thin  curves)  and  0.5  {thick  curves).  Over  the  spec¬ 
tral  index  range  relevant  for  solar  flares,  this  line  ratio  varies  by 
almost  an  order  of  magnitude.  Because  of  the  lower  threshold 
energies  of  a-particle  reactions,  line  production  for  hard  accel¬ 
erated  particle  spectra  are  dominated  by  proton  interactions  and 
there  is  little  difference  between  the  ratios  for  a/proton  =0.1  and 
0.5.  As  the  spectrum  softens,  the  contribution  from  a-particle  in¬ 
teractions  increases  and  the  ratios  differ. 

Accelerated  a-particles  can  interact  with  ambient  He  to  pro¬ 
duce  excited  states  of  7Be  and  7Li.  These  states  then  deexcite  to 
produce  lines  at  0.478  and  0.429  MeV,  which,  through  Doppler 
broadening,  generally  merge  into  what  is  known  as  the  “a-a 
complex.”  Since  the  complex  is  produced  purely  by  accelerated 
a-particles,  the  ratio  of  its  yield  to  the  yield  of  other  deexcita¬ 
tion  lines  (which  are  generally  produced  by  both  protons  and 
a-particles)  is  very  sensitive  to  the  accelerated  a/proton  ratio  in 
addition  to  its  dependence  on  the  accelerated  ion  spectrum.  We 
show  in  Figure  6  the  yields  for  this  complex  as  a  function  of  spec¬ 
tral  index  for  alp  =  0.1  {thin  curves)  and  0.5  {thick  curves).  We 
show  in  Figure  8  the  corresponding  yield  ratios  of  the  a-a  com¬ 
plex  and  the  sum  of  the  strongest  high-FIP  lines  (1.63,4 .44,  6.13, 
and  MeV  lines  from  20Ne,  12C,  160,  and  160,  respectively,  as  a 
function  of  spectral  index.  The  shape  of  this  curve  is  due  to  the  dif¬ 
ferent  energy  dependences  of  the  cross  sections.  The  cross  sections 
for  producing  the  7Li  and  7Be  lines  are  both  sharply  peaked  with 
threshold  energies  of  about  10  MeV  nucleon-1  and  falling  to 
insignificance  by  100  MeV  nucleon-1.  In  contrast,  the  cross  sec¬ 
tions  for  producing  the  4.44  MeV  12C  line  from  both  accelerated 
protons  and  a-particles  have  lower  thresholds  (~5  MeV  for  pro¬ 
ton  reaction  and  1-2  MeV  nucleon-1  for  a  reactions)  and  also  ex¬ 
tend  to  higher  energies,  and  this  is  generally  true  for  the  other  lines 
as  well.  As  a  result,  the  yield  of  the  4.44  MeV  line  is  relatively 
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Fig.  9. — Line  profiles  of  the  4.438  MeV  12C  line  resulting  from  various  in¬ 
teracting-particle  angular  distributions  for  a  disk-centered  flare:  (a)  isotropic, 
( b )  saturated  PAS  and  magnetic  convergence  8  =  0.20,  and  (c)  downward  beam. 
The  solid  and  dotted  curves  are  for  interactions  due  to  protons  and  a-particles, 
respectively.  The  thick  and  thin  curves  are  for  power-law  spectral  indexes  of 
2.75  and  4.75,  respectively.  (Note  that  the  profiles  have  been  renormalized  to 
facilitate  comparison.). 

large  for  steep  spectra  due  to  the  lower  threshold  energy  and  also 
large  for  hard  spectra  due  to  the  extension  to  high  energy,  with  a 
minimum  for  intermediate  spectra. 

2.4. 1.2.  Deexcitation  Line  Shapes  and  Central  Energies 
The  acceleration  parameters  that  most  significantly  affect 
the  shapes  and  line  shifts  of  the  narrow  deexcitation  lines  are  the 
a/proton  ratio  and  the  steepness  of  the  accelerated-particle  spec¬ 
trum  (  the  index  s  for  power-law  spectra).  The  effect  on  the  shape 
of  the  line  is  illustrated  in  Figure  9 a  with  line  profiles  of  the 
4.438  MeV  12C  line  for  an  isotropic  interacting-particle  angular 
distribution.  Because  the  angular  distribution  is  assumed  isotro¬ 
pic,  the  shape  is  symmetrical  about  the  rest  energy  of  the  line 


Fig.  10. — Ratio  of  the  4.438  MeV  12C  line  yields  from  n -particles  and  from 
protons  as  a  function  of  spectral  index  for  equal  numbers  of  a-particles  and 
protons. 

(indicated  by  the  vertical  line).  The  solid  and  dotted  curves  are 
for  interactions  due  to  protons  and  a-particles,  respectively.  The 
thick  and  thin  curves  are  for  power-law  spectral  indexes  of  2.75 
and  4.75,  respectively.  (Note  that  the  profiles  have  been  renor¬ 
malized  to  facilitate  comparison.)  For  both  proton  and  a-particle 
interactions,  a  harder  spectrum  results  in  a  broader  line  because 
there  are  more  interactions  involving  higher  energy  accelerated 
particles  and  such  particles  produce  larger  Doppler  shifts  due  to 
the  faster  recoiling  nuclei.  Similarly,  a-particle  interactions  pro¬ 
duce  a  broader  line  because  the  larger  a-particle  mass  also  re¬ 
sults  in  faster  recoil  velocities. 

The  final  shape  of  the  line  depends  on  the  relative  contribution 
of  the  accelerated  proton  and  a-particle  reactions,  which  depends 
on  the  accelerated  a/proton  ratio  and  how  the  relative  yields  de¬ 
pend  on  spectral  index  through  their  respective  production  cross 
sections.  Because  the  cross  sections  of  a-particle  reactions  tend  to 
have  lower  threshold  energies,  steeper  accelerated-particle  spectra 
tend  to  increase  the  relative  a-reaction  yield.  For  the  4.438  MeV 
12C  line,  Figure  10  shows  the  ratio  of  the  yield  from  a-particles  to 
the  yield  from  protons  as  a  function  of  spectral  index  for  equal 
numbers  of  accelerated  a-particles  and  protons.  Figure  11  shows 
the  final  resulting  line  shapes  for  s  =  3.75  and  a /proton  ratios  of 


Fig.  11. — 4.438  MeV  12C  line  shapes  for  s  =  3.75  and  a/proton  ratios  of 
0.1  ( dotted  line )  and  0.5  {solid  line).  The  accelerated  ion  angular  distribution  is 
isotropic. 
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Fig.  12. — 4.438  MeV  12C  line  profiles  resulting  from  three  angular  distri¬ 
butions  [isotropic,  saturated  PAS  (1  =  20),  and  downward  beam]  for  a  flare  oc¬ 
curring  at  a  heliocentric  angle  of  74.5°  (as  for  the  1991  June  4  flare).  The  o /proton 
ratio  is  0.5  and  s  =  4. 

0.1  and  0.5.  Although  a  harder  spectrum  produces  broader  lines, 
the  contribution  of  the  inherently  broader  a  reaction  decreases 
so  that  the  overall  line  width  may  in  fact  be  narrower  depending 
on  the  a/proton  ratio. 

The  physical  parameters  that  affect  the  shapes  and  shifts  of  the 
narrow  deexcitation  lines  are  those  that  affect  the  angular  distri¬ 
bution  of  the  interacting  particles:  most  significantly  the  level  of 
convergence  (6)  and  PAS  (A)  and  to  a  lesser  degree  the  atmo¬ 
spheric  model.  Figure  9b  illustrates  the  effect  on  the  12C  line  pro¬ 
file  for  6  =  0.20  and  saturated  PAS  (A  =  20)  for  a  flare  centered 
on  the  solar  disk  ( 60bs  =  0°).  Again,  the  solid  and  dashed  curves 
are  for  interactions  due  to  protons  and  a-particles,  respectively, 
and  the  thick  and  thin  curves  are  for  power-law  spectral  indexes 
of  2.75  and  4.75,  respectively.  This  downward-directed  angular 
distribution  (see  Fig.  3)  results  in  redshifted  lines,  with  both  harder 
spectra  and  a-particles  producing  the  most  Doppler  shift.  Details 
of  the  line  shape  and  shift  will  vary  with  different  combinations  of 
6,  A,  s,  and  accelerated  a/proton  ratio.  The  extreme  case  of  a  1 00% 
downward-beamed  angular  distribution  for  a  disk-centered  flare  is 
shown  in  Figure  9c.  As  discussed  in  §  2.3.2,  the  RHESSI/Avretl 
atmosphere  produces  a  slightly  enhanced  downward  component 
for  combinations  of  <5  and  A  resulting  in  mirroring  (i.e.,  fan-beam) 
distributions.  This  results  in  a  slightly  increased  line  redshift  for 
such  distributions.  For  example,  at  the  74.5°  June  4  flare  location 
with  A  =  0  and  S  =  0.45,  the  4.438  MeV  12C  line  is  redshifted 
3  keV  more  for  the  RHESSII Avrett  atmosphere  than  for  the  Avrett 
atmosphere. 

For  any  angular  distribution  with  some  downward-directed 
component,  as  the  flare  location  moves  away  from  Sun  center, 
the  amount  of  Doppler  shift  will  decrease.  This  is  illustrated  in 
Figure  12,  where  12C  line  profiles  are  shown  for  a  flare  occurring 
at  a  heliocentric  angle  of  (90 bs  =  74.5°  (as  for  the  1991  June  4 
flare).  Profiles  resulting  from  three  angular  distributions  are 
shown:  isotropic,  saturated  PAS,  and  downward  beam.  For  these 
calculations,  the  a/proton  ratio  is  0.5  and  s  =  4.  At  74.5°,  only 
the  extreme  case  of  a  100%  downward  beam  produces  any  sig¬ 
nificant  Doppler  shift  of  the  line. 

2.4. 1.3.  Deexcitation  Line  Time  History 

Because  the  lifetimes  of  the  excited  states  are  very  short 
( <  1 0  6  s),  photon  emission  is  prompt,  with  virtually  no  delay 
between  the  nuclear  interaction  and  the  appearance  of  the  7-ray. 
As  a  result,  a  measured  deexcitation  line  emission  time  history 


Fig.  13. — Calculated  yields  of  total  neutrons  and  the  2.223  MeV  neutron- 
capture  line  (for  a  disk-centered  flare)  as  functions  of  the  accelerated-ion  power- 
law  spectral  index  for  a/proton  —  0.5  ( thick  lines)  and  0. 1  ( thin  lines).  The  total 
neutron  yield  has  been  divided  by  a  factor  of  10  for  plotting  purposes.  The 
neutron-capture  curves  are  for  a  disk-centered  flare  and  are  plotted  for  both  no 
PAS  ( dashed  line)  and  saturated  PAS  ( dotted  line).  For  each  PAS,  yields  for 
magnetic  convergences  of  6  =  0.0  and  0.45  are  shown,  as  indicated  in  the  figure. 
All  yields  are  normalized  to  1  accelerated  proton  of  energy  greater  than  30  MeV. 

will  be  essentially  identical  to  the  nuclear  interaction  time  his¬ 
tory.  The  parameters  that  affect  the  interaction  time  history  (and 
thus  the  delay  of  7-ray  emission  relative  to  accelerated  particle 
release)  were  discussed  in  §  2.3.1.  The  only  acceleration  param¬ 
eter  that  affects  the  interaction  time  history  is  the  ion  acceleration 
release  time  history,  a  ;on(t).  The  most  significant  physical  param¬ 
eters  affecting  the  time  history  are  the  magnetic  convergence,  the 
level  of  PAS,  and  the  loop  length. 

2.4.2.  Neutron  Production  and  Escape 

In  this  section  we  discuss  the  measurable  quantities  associated 
with  flare-produced  neutrons  escaping  from  the  Sun  and  observed 
by  detectors  in  space  or  on  Earth.  The  measurable  quantities  are 
the  energies  of  the  neutrons  and  their  arrival  times;  i.e.,  the  time- 
dependent  energy  spectra  of  the  arriving  neutrons.  These  spectra 
are  modified  from  the  spectrum  of  neutrons  that  escaped  from  the 
Sun  by  both  velocity  dispersion  of  the  neutrons  with  various  ener¬ 
gies  and  decay  of  the  neutrons  during  transport.  How  the  accelera¬ 
tion  and  physical  parameters  affect  the  time-dependent  escaping 
neutron  spectrum  was  discussed  by  Hua  et  al.  (2002).  We  sum¬ 
marize  some  of  their  results  here  and  further  discuss  some  other 
aspects  of  neutron  production.  A  fraction  of  those  neutrons  ini¬ 
tially  produced  traveling  in  a  downward  direction  can  be  captured 
on  H  to  produce  deuterium  and  the  2.223  MeV  neutron-capture 
line,  which  is  discussed  in  §  2.4.3. 

2.4.2. 1.  Total  Neutron  Production 

The  acceleration  parameters  that  affect  the  total  neutron  yield 
are  the  accelerated  particle  spectral  index  and  composition,  par¬ 
ticularly  the  accelerated  a/proton  ratio.  The  only  physical  pa¬ 
rameter  affecting  the  total  yield  is  the  ambient  composition.  Total 
neutron  production  as  a  function  of  accelerated-ion  spectral  index 
is  shown  in  Figure  13  for  a/proton  =  0.1  ( thin  curves)  and  0.5 
( thick  curves).  These  total  yields  are  similar  to  those  given  by 
Ramaty  et  al.  (1996)  (who  also  used  the  updated  total  cross  sec¬ 
tions  of  Hua  et  al.  2002),  except  for  steep  spectra  where  the  yields 
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TABLE  3 

Approximate  Threshold  Energies  of  Neutron 
Production  Reactions 


Reaction 

Threshold  Energy 
(MeV  nucleon-1) 

p-CNO . 

5-10 

cif-CNO . 

~1 

V  'I . 

~300 

/>  lie . 

~30 

a-He . 

~10 

calculated  here  are  ~30%  larger  due  to  our  larger  assumed  am¬ 
bient  He/H  abundance.  Using  the  loop  formalism  of  Hua  et  al. 
(1989, 2002)  we  confirmed  that  this  total  neutron  yield  is  indepen¬ 
dent  ofparameters  of  the  loop  model  as  expected.  For  softer  spec¬ 
tra,  neutron  production  is  dominated  by  a-particle  reactions  and 
the  dependence  of  the  yield  on  the  a/proton  ratio  is  almost  linear. 
The  explanation  for  the  shapes  of  the  yield  curves  is  the  same  as 
for  the  deexcitation  lines,  which  is  discussed  in  §  2.4.1. 

2.4.2.2.  Neutron  Spectra 

In  the  analysis  of  neutrons  from  solar  flares,  there  are  three 
relevant  energy  spectra  of  neutrons  that  must  be  considered:  the 
spectrum  at  the  production  site,  the  spectrum  escaping  from  the 
Sun,  and  the  spectrum  arriving  at  Earth.  As  for  total  neutron  pro¬ 
duction,  the  angle-integrated  spectrum  of  neutrons  at  the  produc¬ 
tion  site  is  independent  of  parameters  of  the  loop  model,  depending 
only  on  the  accelerated-particle  spectrum  and  the  ambient  and 
accelerated-ion  compositions.  Neutrons  are  produced  in  proton- 
proton,  a-a,  proton-a,  and  a-proton  reactions  and  reactions  of 
protons  and  a-particles  with  heavier  ambient  nuclei  and  their  in¬ 
verse  reactions.  The  threshold  energies  for  these  various  reac¬ 
tions  (see  Hua  et  al.  (2002)  for  the  full  cross  sections)  are  given  in 
Table  3  and  span  a  broad  range  of  energies  from  about  an  MeV 
nucleon-1  for  reactions  involving  heavy  nuclei  to  several  hun¬ 
dred  MeV  for  proton-proton  reactions.  Which  reactions  domi¬ 
nate  neutron  production  at  a  given  neutron  energy  depends  on 
the  accelerated-ion  spectrum  and  the  composition  of  the  acceler¬ 
ated  ions  and  the  ambient  medium. 

Angle-integrated  spectra  at  the  production  site  are  shown  in 
Figure  14,  calculated  for  an  accelerated-ion  power- law  spectral 
index  of  4,  saturated  PAS,  6=  0.2,  the  compositions  of  Table  1, 
and  both  accelerated  a/proton  =  0.5  (Fig.  14a)and0.1  (Fig.  14  b). 
Also  shown  are  the  contributions  from  the  various  reactions.  The 
shapes  of  the  total  spectra  for  the  two  values  of  a/proton  are  sim¬ 
ilar,  but  the  overall  yield  is  less  for  a/proton  =  0.1,  as  shown  in 
Figure  13.  For  this  spectral  index,  the  main  contribution  to  the 
production  of  high-energy  (>1  GeV  )  neutrons  is  from  acceler¬ 
ated  a-particle  reactions  with  ambient  H  for  a/proton  =  0.5,  but 
it  is  from  accelerated  proton  reactions  with  ambient  H  and  He  for 
a/proton  =  0.1.  For  both  values  of  a/proton,  the  main  contribu¬ 
tion  to  production  of  neutrons  with  energies  of  tens  of  MeV  is 
from  a-particle  reactions  with  ambient  H.  Referring  to  the  thresh¬ 
olds  of  Table  2,  these  a-particles  must  have  energies  greater  than 
10  MeV  nucleon- '.For  accelerated  a/proton  =  0.5,  the  main  con¬ 
tribution  to  low-energy  (<  10  MeV  )  neutrons  is  from  a  reactions 
with  ambient  He  and  heavier  nuclei.  Because  of  the  very  low 
threshold  for  the  latter  reactions,  some  of  these  low-energy  neu¬ 
trons  can  be  produced  by  a-particles  of  only  a  few  MeV  nucleon- 1 . 
For  a/proton  =0.1,  most  of  the  low-energy  neutrons  are  pro¬ 
duced  by  accelerated  proton  interactions  with  ambient  He.  From 
Table  2,  the  energies  of  these  protons  must  be  >30  MeV.  As  the 
assumed  accelerated-ion  spectrum  steepens,  the  contribution  of 


Fig.  14. — Angle-integrated  neutron  energy  spectra  calculated  for  an  accelerated- 
ion  power-law  spectral  index  s  —  4,  saturated  PAS,  6  =  0.2,  the  compositions 
given  in  Table  1,  and  accelerated  a  /proton  =  0.5  (a)  and  0.1  ( b ).  Contributions 
to  the  total  production  by  each  of  the  various  reactions  are  shown.  All  yields  are 
normalized  to  1  accelerated  proton  of  energy  greater  than  30  MeV. 


the  low-threshold  a-heavy  reactions  begins  to  dominate.  We  note 
that  the  combined  yield  of  neutrons  from  reactions  of  accelerated 
protons  with  ambient  C  and  heavier  nuclei  and  their  inverse  re¬ 
actions  (Fig.  14,  dash-dot-dotted  curves)  increases  as  accelerated 
a/proton  is  increased  even  though  the  reactions  do  not  involve 
accelerated  a-particles.  This  is  because  we  maintain  accelerated 
a/O  at  50,  as  stated  in  §  2.2.1,  which  affects  the  yields  of  the  in¬ 
verse  reactions.  The  relative  contribution  of  the  various  reactions 
to  total  neutron  production  is  independent  of  the  loop  parameters 
such  as  PAS,  magnetic  convergence,  and  length. 

Because  most  of  the  neutrons  are  produced  at  densities  where 
scattering  is  small  (see  Fig.  9  of  Hua  et  al.  2002),  neutrons  ini¬ 
tially  produced  moving  upward  typically  escape  from  the  solar 
atmosphere  without  scattering.  Neutrons  moving  downward  typ¬ 
ically  lose  their  energy  with  little  chance  of  escape.  As  discussed 
by  Hua  et  al.  (2002),  while  low-energy  neutrons  tend  to  be  pro¬ 
duced  isotropically,  higher  energy  neutrons  tend  to  be  emitted 
more  in  the  same  direction  as  the  interacting  ion,  with  this  ten¬ 
dency  increasing  with  increasing  ion  energy.  As  a  result,  the  spec¬ 
trum  of  neutrons  escaping  from  the  Sun  along  a  particular  line  of 
sight  depends  on  the  energy  and  angular  distributions  of  the  in¬ 
teracting  accelerated  ions.  As  discussed  in  §  2.3.2,  the  physical 
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Fig.  15. — Angle-dependent  escaping  neutron  spectra  at  the  Sun.  For  loops 
perpendicular  to  the  solar  surface,  9  is  the  angle  relative  to  the  surface  normal, 
with  6  —  0°  being  the  direction  away  from  the  Sun,  and  corresponds  to  the 
heliocentric  angle  of  a  flare,  0Obs-  For  a  disk-centered  flare,  cos  (9)  =  1.  (a)  No 
PAS  ( solid  curves )  and  saturated  PAS  ( dashed  curves )  with  6  =  0.20.  ( b )  Mag¬ 
netic  convergence  8  =  0.20  ( solid  curves )  and  0.45  ( dashed  curves )  with  A  = 
300.  All  yields  are  normalized  to  1  accelerated  proton  of  energy  greater  than 
30  MeV. 


parameters  that  affect  the  angular  distribution  of  the  interacting 
ions  are  the  level  of  convergence  (8)  and  PAS  (A).  In  Figure  15a 
we  show  the  angular-dependent  kinetic-energy  spectra  of  escap¬ 
ing  neutrons  calculated  for  an  accelerated-ion  spectral-index  of  4, 
convergence  8  =  0.2,  and  both  saturated  (A  =  20,  dashed  line) 
and  no  (A  — >  oo,  solid  line)  PAS.  For  loops  perpendicular  to  the 
solar  surface,  9  is  the  angle  between  the  normal  to  the  solar  surface 
and  the  line  of  sight  and  corresponds  to  the  heliocentric  angle  of 
the  flare,  90 bs.  At  each  angle,  the  number  of  escaping  neutrons 
from  saturated  PAS  is  less  than  that  from  no  PAS  since  fewer  neu¬ 
trons  are  produced  moving  upward  with  saturated  PAS. 

In  Figure  15 b  we  show  the  angular-dependent  kinetic-energy 
spectra  of  escaping  neutrons  calculated  for  an  accelerated  ion 
spectral  index  of  4,  A  =  300,  and  two  values  of  magnetic  conver¬ 
gence,  8  =  0.20  ( solid  line)  and  0.45  ( dashed  line).  The  effect  is 
similar  to  that  for  PAS,  with  the  more  downward-directed  angu¬ 
lar  distribution  of  8  =  0.2  resulting  in  fewer  neutrons  directed 
upward.  The  broader  interacting-ion  angular  distribution  of  8  = 
0.45  (see  Fig.  3)  results  in  a  weaker  dependence  of  the  escaping 
neutron  spectrum  on  angle. 

2.4. 2. 3.  Neutron  Production  Time  History 

The  effects  of  the  acceleration  and  physical  parameters  on  the 
delay  of  the  production  of  neutrons  relative  to  the  accelerated 
ion  release  are  qualitatively  the  same  as  those  described  for  the 
4.44  MeV  12C  deexcitation  line  (§  2.2.1).  The  effects  of  PAS  on 


the  time  dependence  of  neutron  production  were  specifically  dis¬ 
cussed  by  Hua  et  al.  (2002). 

2.4.3.  Neutron-Capture  Line  Production 

In  this  section  we  discuss  the  measurable  quantities  associated 
with  the  2.223  MeV  neutron  capture  line:  the  line  yield  and  time 
history.  Neutrons  initially  produced  moving  downward  either 
decay,  react  with  3  He,  or  are  captured  on  H.  Capture  on  other  nu¬ 
clei  is  less  important  due  to  their  smaller  relative  abundances. 
The  reaction  with  3  He  is  charge  exchange,  3He(«,/>)3H,  without 
the  emission  of  radiation.  The  capture  on  H  results  in  the  formation 
of  deuterium  with  the  binding  energy  appearing  as  a  2.223  MeV 
photon.  Since  the  probability  for  elastic  scattering  is  much  larger 
than  the  probability  for  either  of  these  reactions,  most  of  the 
neutrons  thermalize  first,  causing  a  delay  in  the  formation  of  the 
capture  photon  relative  to  neutron  formation.  This  delay  is  in  ad¬ 
dition  to  the  delay  of  neutron  production  relative  to  the  accel¬ 
erated  ion  release  discussed  in  §  2.4.2. 

The  cross  sections  for  the  H  and  3He  reactions  are  2.2  x 
nr30/?-1  cm2  and  3.7  x  ICC26/?-1  cm2,  respectively,  where  (3  is 
neutron  velocity  in  units  of  the  speed  of  light.  Therefore,  if  the 
3He/H  ratio  at  the  capture  site  is  ~5  x  lCV5  (comparable  to  that 
observed  in  the  solar  wind),  nearly  equal  numbers  of  neutrons 
will  be  captured  on  H  as  react  with  3He.  The  ambient  3He/ H  abun¬ 
dance  ratio  affects  not  only  the  total  yield  of  the  capture  line  but 
also  the  delay  of  its  formation. 

In  the  next  section  we  calculate  2.223  MeV  neutron-capture 
line  yields  and  in  the  following  section  neutron-capture  line  time 
histories.  We  explore  their  dependences  on  the  acceleration  and 
physical  parameters. 

2.4.3. 1.  Neutron-Capture  Line  Yield 

The  total  neutron  capture  line  yield  depends  directly  on  the  to¬ 
tal  neutron  yield  and  so,  in  addition  to  ambient  3He/H,  it  also  de¬ 
pends  on  the  same  acceleration  and  physical  parameters  affecting 
total  neutron  production:  accelerated  particle  spectral  index  and 
ambient  and  accelerated  particle  compositions.  Because  the  prob¬ 
ability  for  capture  depends  on  the  initial  direction  of  the  neutron 
through  the  atmosphere,  the  capture  line  yield  also  depends  on  the 
physical  parameters  affecting  the  interacting  ion  angular  distribu¬ 
tion,  8  and  A,  because  it  determines  the  subsequent  neutron  an¬ 
gular  distribution.  Because  effective  capture  requires  high  density, 
the  capture  line  is  produced  deep  in  the  photosphere,  resulting  in 
strong  attenuation  of  the  line,  especially  if  the  flare  is  located  near 
the  solar  limb.  Therefore,  while  deexcitation  line  yields  do  not  de¬ 
pend  on  the  location  of  the  flare  on  the  solar  surface  (see  §2.4.1), 
the  neutron-capture  line  yield  does. 

Capture-line  yields  for  a  flare  occurring  near  the  center  of  the 
solar  disk  are  shown  in  Figure  13  as  functions  of  the  power-law 
spectral  index  for  a/proton  =0.1  (thin  curves)  and  0.5  ( thick 
curves).  For  each  value  of  a/proton,  the  yields  are  shown  for  two 
extreme  values  of  PAS:  no  PAS  (A  — >  oo,  dashed  curves)  and 
saturated  PAS  (A  =  20,  dotted  curves).  For  each  of  these  PAS 
values,  yields  for  two  extreme  values  of  magnetic  convergence 
are  shown:  8  =  0.0  (i.e.,  none)  and  8  =  0.45.  For  these  calcula¬ 
tions,  the  loop  length  is  11,500  km,  the  atmosphere  is  Avrett 
(1981)and3He/H  =  3.7  x  1 0  5.  As  for  the  deexcitation  line  and 
neutron  yields,  the  shape  of  a  curve  results  from  (1)  the  relation¬ 
ship  between  the  spectrum  and  the  energy  dependence  of  the  cross 
section  and  (2)  the  normalization  of  1  proton  above  30  MeV.  We 
see  that  the  yield  is  larger  for  saturated  PAS  than  for  no  PAS  since 
the  more  downward-directed  neutrons  from  PAS  are  more  likely 
to  be  captured  than  decay.  Similarly,  the  yields  are  larger  for  smaller 
values  of  convergence,  again  due  to  the  more  downward-directed 
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angular  distribution.  We  note  that  this  is  in  contrast  to  the  es¬ 
caping  neutron  yield  discussed  in  §  2.4.2,  which  is  less  for  more 
downward-directed  angular  distributions. 

The  ambient  3He/H  ratio  has  a  strong  effect  on  the  neutron- 
capture  line  yield  as  the  nonradiative  reactions  with  'He  compete 
with  the  line-producing  captures  on  H.  As  3He/H  is  increased 
from  3.7  x  10~5  to  10  x  10  5,  the  yield  decreases  by  ^50%  and 
when  3He/H  is  decreased  to  0,  the  yield  increases  by  ~50%. 
Up  to  heliocentric  angles  of  85°,  the  dependence  of  the  yield  on 
3He/H  is  fairly  independent  of  interacting  particle  angular  distri¬ 
bution  (i.e.,  X  and  8),  spectral  index,  accelerated  a/proton  ratio, 
and  atmospheric  model. 

We  found  that  the  effect  of  the  atmospheric  model  on  the  total 
yield  was  relatively  weak,  with  the  most  significant  differences 
occurring  for  combinations  of  8  and  X  producing  mirroring  (i.e., 
fan-beam)  distributions  of  the  interacting  ions  (and  resulting  neu¬ 
trons).  The  fate  of  downward-moving  neutrons  is  fairly  insensi¬ 
tive  to  the  atmospheric  structure,  but  the  relative  probabilities  for 
capture  and  decay  for  neutrons  moving  parallel  to  the  surface  de¬ 
pend  on  the  amount  of  material  the  neutrons  encounter.  The  higher 
density  of  the  RHESSII Avrett  model  at  the  higher  elevations 
resulted  in  higher  capture  line  yields  for  the  fan-beam  angular 
distributions.  For  parameter  combinations  resulting  in  strongly 
downward  distributions,  the  yields  from  the  Avrett  and  RHESSI/ 
Avrett  models  were  the  same,  regardless  of  flare  location.  The 
yields  from  the  Vemazza  model  were  as  much  as  1 0%  less,  with 
the  largest  difference  occurring  for  limb  flares.  For  parameter 
combinations  resulting  in  more  fan-beam  distributions,  the  yields 
from  the  Avrett  model  were  ~  1 0%  less  than  from  the  RHESSI 7 
Avrett  model,  independent  of  flare  location.  The  yields  from  the 
Vemazza  model  were  from  15%  to  20%  less  than  the  RHESSI/ 
Avrett  model,  independent  of  heliocentric  angle.  These  effects 
were  not  significantly  dependent  on  3He/H. 

Because  of  the  very  different  energy  dependences  of  deexcita¬ 
tion  line  and  neutron  production  cross  sections,  the  fluence  ratio 
of  the  neutron-capture  line  and  a  deexcitation  line  is  very  sensi¬ 
tive  to  the  accelerated-ion  spectral  index.  Corresponding  to  the 
yields  shown  in  Figure  13,  neutron-capture/4.438  MeV  12C  line 
fluence  ratios  as  a  function  of  power-law  spectral  index  for  a  flare 
occurring  near  the  center  of  the  solar  disk  are  shown  in  Figure  7 
for  a/proton  =0.1  ( thin  curves )  and  0.5  ( thick  curves).  For  each 
value  of  a/proton,  the  yields  are  shown  for  two  extreme  values  of 
PAS:  no  PAS  ( X  — >  oo,  dashed  curves)  and  saturated  PAS  ( X  = 
20,  dotted  curves).  For  each  of  these  PAS  values,  yields  for  two 
extreme  values  of  magnetic  convergence  are  shown:  6  —  0.0 
(i.e.,  none)  and  6  =  0.45.  The  ratios  vary  by  more  than  2  orders 
of  magnitude  as  the  index  is  varied  and  are  significantly  more 
sensitive  and  cover  a  wider  accelerated-ion  energy  range  than  the 
6.129  MeV  160/1.634  MeV  20Ne  ratio  (see  §  2.4.4).  Since  the 
neutron-capture  line  is  delayed  by  ~  1 00  s  due  to  the  time  required 
by  the  neutrons  to  slow  down  and  be  captured  (see  below),  this 
ratio  generally  is  only  used  to  provide  flare-averaged  spectral  in¬ 
dexes.  For  these  calculations,  3He/H  =  3.7  x  1 0  5 . 

The  neutron-capture  line  originates  in  the  photosphere  (in 
contrast  to  nuclear  deexcitation  lines,  which  probably  origi¬ 
nate  at  chromospheric  densities)  and  is  significantly  attenuated 
for  flares  that  occur  near  the  solar  limb.  We  show  in  Figure  16 
neutron-capture/4.438  MeV  12C  line  ratios  as  functions  of  flare 
heliocentric  angle,  parameterized  by  spectral  index,  and  calcu¬ 
lated  for  accelerated  a/proton  =  0.5, 3He/H  =  3  x  10-5,  and  no 
PAS  (X  — >  oo;Fig.  1 6a),  nearly  saturated  PAS  (1  =  300;  Fig.  16 b), 
and  saturated  PAS  (X  =  20;  Fig.  16c).  The  neutron-capture  line 
is  strongly  attenuated  for  flares  occurring  at  heliocentric  angles 
greater  than  ^75°:  relative  to  a  disk-centered  flare,  the  neutron- 


Fig.  16. — The  2.223  MeV  (neutron  capture)/4.438  MeV  ( 12C)  line  fluence  ratio 
as  a  function  of  the  flare  heliocentric  angle  (0Obs-)>  parameterized  y  spectral  index 
(5):  (a)  no  PAS,  (b)  nearly  saturated  PAS,  and  (c)  saturated  PAS. 

capture/4.438  MeV  12  C  line  ratio  decreases  by  an  order  of  mag¬ 
nitude  by  88°  and  2  orders  of  magnitude  by  90°.  We  estimate  the 
uncertainty  of  the  neutron-capture  yield  due  to  uncertainty  in  the 
neutron  production  cross  sections  to  be  ~20%. 

2.4. 3. 2.  Neutron-Capture  Line  Time  History 

Because  both  the  probability  for  capture  and  the  time  to  ther- 
malize  before  capture  depend  on  the  initial  direction  of  the  neu¬ 
tron  through  the  atmosphere,  the  capture  line  time  history  depends 
on  the  physical  parameters  affecting  the  interacting  particle  an¬ 
gular  distribution:  X  and  8.  In  Figure  1  la  we  show  neutron-capture 
line  time  histories  calculated  for  instantaneous  release  of  acceler¬ 
ated  ions  at  the  top  of  the  loop  at  t  =  0  for  various  values  of  PAS, 
from  X  — >  00  (none)  to  X  =  20  (saturated).  (We  note  that  these 
time  histories,  and  the  other  time  histories  shown  in  this  section, 
include  the  delay  of  neutron  production  relative  to  the  initial  ion 
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Fig.  17. — Dependence  of  the  2.223  MeV  neutron-capture  line  time  history  on  various  parameters:  (a)  PAS  (2);  ( b )  magnetic  convergence  (6);  (c)  spectral  index 
(5);  ( d )  ambient  3He/H;  ( e )  flare  heliocentric  angle  (0Obs);  and  (/)  atmospheric  model.  For  each  panel,  the  other  parameters  are  fixed  at  the  values  noted  in  the  text.  All 
yields  are  normalized  to  1  accelerated  proton  of  energy  greater  than  30  MeV. 


release.)  For  these  calculations,  accelerated  a/proton  =  0.5, 
3He/H  =  3.7  x  10-5,  and  0obs  =  75°  (where  0o bS  is  the  heliocen¬ 
tric  angle  of  the  flare).  As  PAS  decreases  (A  increases),  the 
interacting-particle  angular  distribution  (and  corresponding  neu¬ 
tron  distribution)  become  less  downward-directed  and  the  time 
history  becomes  more  prolonged.  Similarly,  the  angular  distribu¬ 
tion  becomes  less  downward  directed  as  6  decreases  and  the  time 
history  again  becomes  more  prolonged,  as  shown  in  Figure  17 b 
(for  A  =  300). 

Because  harder  accelerated  particle  spectra  produce  more  high- 
energy  neutrons,  which  penetrate  to  higher  densities  with  in¬ 
creased  probability  for  capture,  the  capture  line  time  history  also 
depends  on  the  accelerated  particle  spectral  index.  This  is  shown 
in  Figure  1 7c  for  several  values  of  spectral  index  s.  Harder  spectra 
result  in  a  faster  falling  time  history.  For  these  calculations,  ac¬ 
celerated  a/proton  =  0.5, 3He/H  =  3.7  x  10~5,  and  0obs  =  75°. 
Also  apparent  is  the  strong  dependence  of  the  total  line  yield  on 
spectral  index. 

The  nonradiative  reaction  with  3He  competes  with  neutron 
capture  on  H  and  affects  the  neutron-capture  line  time  history. 
This  is  seen  in  Figure  1  Id,  where  time  histories  for  several  values 


of  3He/H  are  shown.  Increased  3He  results  in  a  faster  falling  time 
history  as  neutrons  are  more  rapidly  removed.  For  these  calcu¬ 
lations,  accelerated  a/proton  =  0.5  and  0obs  =  75°. 

Neutrons  that  penetrate  deeper  in  the  photosphere  are  captured 
sooner  due  to  the  higher  density,  but  the  subsequent  capture  line 
photons  are  more  strongly  attenuated.  For  flares  located  nearer  to 
the  solar  limb,  this  attenuation  becomes  even  stronger,  resulting 
in  slower  time  histories  as  photons  from  the  more-rapid  captures 
are  attenuated.  This  can  be  seen  in  Figure  17e,  where  we  show 
capture-line  time  histories  for  flares  located  at  disk  center  (0°  helio¬ 
centric  angle),  15° ,  and  85°.  Also  apparent  is  the  strong  decrease 
of  the  total  line  yield  with  increasing  angle.  For  these  calcula¬ 
tions,  accelerated  a/proton  =  0.5  and  3He/H  =  3.7  x  1CU5. 

We  found  that  the  effect  of  the  atmospheric  model  on  the 
neutron-capture  line  time  history  was  relatively  weak,  with  the 
most  significant  differences  occurring  for  combinations  of  6  and 
A  resulting  in  interacting  particle  (and  resulting  neutron)  angular 
distributions  that  were  the  least  downward  directed.  Similar  to 
the  total  capture  line  yield  as  noted  above,  the  fate  of  downward- 
moving  neutrons  is  fairly  insensitive  to  the  atmospheric  structure, 
but  the  relative  probabilities  for  capture  and  decay  for  neutrons 
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Fig.  18. — 1.634  MeV  20Ne  deexcitation  line,  2.223  MeV  neutron  capture 
line  and  positron  yields  per  accelerated  proton  of  given  energy.  For  the  Ne  line, 
the  yield  is  total  photons,  for  the  neutron  capture  line,  the  yield  is  photons  sr_  1  at 
heliocentric  angle  f?0bs  =  75°,  and  for  positrons,  the  yield  is  total  positrons. 

moving  parallel  to  the  surface  depend  on  the  amount  of  material 
the  neutrons  encounter.  The  higher  density  of  the  RHESSI/ Avrett 
model  at  the  higher  elevations  results  in  a  faster  decay  of  the  time 
history  for  fan-beam  distributions.  This  can  be  seen  in  Figure  17 'f 
where  we  show  time  histories  calculated  for  the  Avrett  and  the 
RHESSI/Avrett  atmospheres.  Here  we  assume  6  =  0.2  and  A  = 
1 000,  which  results  in  a  fan  beam  angular  distribution  (see  Fig.  3c). 
For  the  other  parameters,  accelerated  a/proton  =  0.5,3He/H  = 
3.7  x  10-5,  and  0obs  =  75°. 

Where  neutrons  are  formed  in  the  loop  does  not  depend  on 
loop  length  (see  §  2.3.3),  so  the  time  required  forneutrons  to  sub¬ 
sequently  thermalize  and  be  captured  also  does  not  depend  on 
loop  length.  The  only  effect  of  loop  length  on  the  neutron-capture 
line  time  history  is  due  to  the  time  dependence  of  the  neutron 
production  itself  (see  §  2.3 . 1 ).  We  find  that,  unless  the  loop  is  quite 
long  (~80,000  km),  the  effect  is  limited  to  only  the  first  1 0—30  s, 
producing  a  slight  delay  of  the  neutron-capture  time  history  in¬ 
creasing  to  ~30  s  with  increasing  loop  length.  After  about  30  s, 
the  capture  line  time  history  is  unaffected. 

2.4.4.  Effective  Energies  of  the  Accelerated  Ions 
Producing  Gamma  Rays  and  Neutrons 

In  the  previous  sections  we  have  explored  the  dependences  of 
the  various  measurable  quantities  of  flare  emissions  on  the  accel¬ 
eration  and  physical  parameters.  When  deriving  these  parame¬ 
ters  from  measured  observables,  it  is  important  to  be  aware  that 
each  observable  is  sensitive  only  to  a  unique  range  of  accelerated 
ion  energies.  For  example,  when  deriving  spectral  indexes  from 
measured  line  fluence  ratios,  the  effective  ion  energies  relevant 
to  the  16O/20Ne  ratio  are  much  lower  than  those  relevant  to  the 
neutron  capture/12C  line  ratio.  In  this  section  we  determine  the 
effective  ion  energy  ranges  responsible  for  the  various  products 
of  nuclear  reactions. 

The  main  factor  determining  the  effective  energy  range  for  a 
nuclear  reaction  is  the  cross  section,  but  the  shape  of  the  acceler¬ 
ated  ion  spectrum  can  be  important.  Steep  spectra  can  shift  the 
emphasis  to  lower  energies  and  vice  versa.  For  thin  target  inter¬ 
actions,  only  the  cross  section  and  spectrum  are  relevant.  As  a 
general  rule,  only  ions  with  energies  where  the  cross  section  is 
significant  contribute  to  thin  target  interactions.  For  thick  target 
interactions,  however,  even  ions  with  high  energies  where  cross 
sections  are  often  small  also  contribute,  since  they  can  lose  energy 
and  then  interact  where  the  cross  section  is  significant. 


Proton  energy  (MeV) 

Fig.  19. — Yield  of  the  6.129  MeV  160  line  weighted  by  accelerated  ion 
power-law  spectra  with  indexes  5  =  3  and  5.  The  contribution  to  the  line  from 
accelerated  protons  and  a-particles  are  separately  shown. 

Using  our  nuclear  production  computer  codes  we  have  cal¬ 
culated  the  yields  of  deexcitation  lines,  escaping  neutron  spec¬ 
tra  and  the  neutron  capture  line  for  monoenergetic  accelerated 
ions.  For  the  deexcitation  lines  we  include  all  spallation  reactions 
leading  to  the  excited  nuclei.  We  show  in  Figure  18  monoener¬ 
getic  ion  yields  for  the  1 .634  MeV  20Ne  line  and  the  2.223  MeV 
neutron  capture  line.  The  composition  of  the  accelerated  particles 
is  that  given  in  Table  2,  and  the  yields  are  normalized  to  one  pro¬ 
ton  at  the  given  energy.  Accelerated  ions  heavier  than  He  contrib¬ 
ute  to  the  neutron  capture  line,  but  for  the  deexcitation  lines  we  are 
only  considering  accelerated  proton  and  a  particle  direct  reactions 
producing  the  narrow  lines.  For  the  neutron  capture  line  calculation, 
we  assumed  A  =  300,  6  =  0.2,  3He/H  =  3  x  10  5,  and  0o bS  = 
75°.  The  effect  of  the  depth  distribution  of  the  production  and 
subsequent  attenuation  of  the  2.223  MeV  photons  has  been  taken 
into  account.  Because  of  the  thick  target  assumption,  we  see  that 
ions  of  several  hundred  MeV  nucleon- 1  are  very  effective  in  pro¬ 
ducing  the  Ne  line  even  though  the  cross  section  at  those  energies 
is  small  (see  Kozlovsky  et  al.  2002). 

While  we  do  not  discuss  in  detail  the  production  of  positrons 
in  this  paper,  the  positron  annihilation  line  is  an  important  addi¬ 
tional  diagnostic  for  flare  parameters.  We  also  show  in  Figure  1 8 
positron  yields  for  monoenergetic  ions  and  emphasize  that  these 
yields  are  for  positrons,  not  the  subsequent  annihilation  photons. 
The  assumed  accelerated  3He/4He  ratio  is  1.0.  The  structure  ap¬ 
parent  in  this  curve  deserves  discussion.  The  enhancement  at 
high  ion  energies  (>300  MeV  nucleon-1)  is  due  to  positrons  from 
the  decay  of  positively  charged  pions  whose  production  cross  sec¬ 
tions  have  very  high  threshold  energies  (Murphy  et  al.  1987).  The 
enhancement  at  low  energies  (less  than  ^5  MeV  nucleon-1)  is 
due  to  positrons  from  the  decay  of  radioactive  nuclei  produced  by 
accelerated  3He  reactions,  which  can  have  large  cross  sections  at 
very  low  energies  (Kozlovsky  et  al.  2004).  The  yield  at  interme¬ 
diate  energies  is  due  to  positrons  from  the  decay  of  radioactive  nu¬ 
clei  produced  by  accelerated  protons  and  a-particles  with  heavier 
elements  and  their  inverse  reactions  (Kozlovsky  et  al.  1987). 

The  differential  accelerated  ion  kinetic  energy  spectrum  gives 
the  number  of  ions  per  unit  energy  at  each  energy.  To  deter¬ 
mine  the  effective  ion  energies  for  a  particular  ion  spectrum,  we 
weight  the  monoenergetic  yields  calculated  above  with  the  given 
spectrum.  As  an  example,  we  show  in  Figure  19  the  differen¬ 
tial  yield  of  the  O  line  as  a  function  of  accelerated  ion  energy 
for  power-law  spectra  with  indexes  of  3  and  5  and  accelerated 
a/proton  =  0.5.  We  show  the  separate  contributions  to  the  line 


184 


MURPHY  ET  AL. 


Vol.  168 


Spectral  index 

Fig.  20. — Effective  accelerated  ion  energies  for  production  of  the  1 .634  MeV 
20Ne  line  ( white  boxes)  and  6.129  MeV  160  line  ( gray  boxes)  as  a  function  of 
spectral  index.  The  effective  energy  ranges  are  the  50%  yield  range  as  defined  in 
the  text,  (a)  Accelerated  a/proton  =  0.1;  (b)  accelerated  a/proton  =  0.5. 

from  accelerated  protons  ( dashed  curves )  and  a-particles  ( dotted 
curves).  We  see  that  for  a  soft  spectrum,  the  a  particle  contribu¬ 
tions  dominate  due  their  lower  threshold  energies,  and  the  most 
effective  ion  energies  are  around  a  few  MeV  nucleon-1.  For  hard 
spectra,  the  proton  interactions  become  important  and  the  effec¬ 
tive  ion  energies  shift  to  higher  energies  due  to  the  higher  proton 
interaction  threshold  energies.  For  intermediate  spectral  indexes, 
both  accelerated  particles  can  be  important  and  the  effective  ion 
energy  range  can  be  very  broad.  To  quantify  the  effective  ion  en¬ 
ergy  range  for  producing  the  line,  we  use  the  following  definition: 
we  identify  the  ion  energy  where  the  yield  is  maximum  and  then 
define  the  effective  range  to  be  that  within  which  the  yield  has 
fallen  to  50%  on  each  side  of  the  maximum.  We  note  that  if  the 
effective  energy  distribution  is  very  broad,  the  arbitrary  value  of 
50%  could  be  somewhat  misleading. 

The  fluence  ratio  of  the  6.129  MeV  160  and  1.634  MeV  20Ne 
lines  has  been  frequently  used  as  a  measure  of  the  accelerated  ion 
spectral  index.  Using  the  50%  definition  for  the  effective  accel¬ 
erated  ion  energy  range  as  given  above,  we  can  now  see  what  ion 
energies  are  most  relevant  to  this  ratio.  We  show  in  Figure  20 
the  effective  accelerated  ion  energy  ranges  for  producing  the 
1 .634  MeV  20Ne  line  ( white  boxes)  and  the  6. 129  MeV  160  line 
( gray  boxes)  for  several  power-law  spectral  indexes  and  for  ac¬ 
celerated  ct/proton  =  0.1  (Fig.  20a)  and  0.5  (Fig.  20b).  The 
horizontal  line  within  each  box  corresponds  to  the  peak  of  the 
distribution.  The  ratio  of  these  two  lines  has  been  used  to  de¬ 
termine  accelerated  ion  spectral  indexes  (see  §  2.4.1).  The  fig¬ 
ure  shows  that  the  effective  energy  range  for  producing  the  160 
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Fig.  2 1 . — Effective  accelerated  ion  energies  for  production  of  the  4.43  8  MeV 
12C  line  ( white  boxes)  and  2.223  MeV  neutron  capture  line  ( gray  boxes)  as  a 
function  of  spectral  index.  The  effective  energy  ranges  are  the  50%  yield  range 
as  defined  in  the  text.  Panel  (a)  shows  accelerated  a/proton  =  0.1  and  (b)  shows 
a/proton  =  0.5.  For  the  neutron  capture  line,  X  =  300,  8  =  0.2,  L  =  1 1,500  km, 
and  6Us  =  75°. 

line  is  always  shifted  to  higher  energies  than  for  the  20Ne  line. 
This  is  reflected  in  Figure  7,  where  the  16O/20Ne  yield  ratio 
curve  monotonically  decreases  with  increasing  spectral  index. 
The  degree  of  separation  of  the  effective  energy  ranges  for  the 
two  lines  determines  the  sensitivity  of  the  ratio  to  the  spectral 
index  (i.e.,  the  slope  of  the  curve  in  Fig.  7).  Because  the  separa¬ 
tion  for  these  two  lines  is  not  large,  the  slope  of  the  yield  ratio 
curve  in  Figure  7  is  not  very  steep. 

The  range  of  accelerated  ion  energies  relevant  to  this  line  ratio 
can  also  be  seen  in  Figure  20.  For  very  hard  spectra,  the  upper 
range  of  the  effective  ion  energies  can  be  ~  100  MeV  nucleon-1. 
But  for  most  (softer)  spectra  it  is  less  than  10  MeV  nucleon-1. 
The  lower  range  is  typically  a  few  MeV  nucleon-1,  but  for  very 
soft  spectra  (particularly  when  a/proton  =  0.5)  it  can  be  less  than 
2  MeV  nucleon-1.  For  softer  spectra,  the  line  production  is  dom¬ 
inated  by  the  lower  energy  reactions  of  accelerated  a-particles. 
As  the  spectrum  hardens,  the  higher  energy  accelerated  proton 
reactions  begin  to  contribute  and  the  effective  ion  energy  range 
shifts  to  higher  energies  and  becomes  broader  as  interactions  of 
both  a-particlcs  and  protons  contribute. 

The  fluence  ratio  of  the  4.438  MeV  12C  and  2.223  MeV  neu¬ 
tron  capture  lines  has  also  been  frequently  used  as  a  measure  of 
the  accelerated  ion  spectral  index.  We  show  in  Figure  21  the  ef¬ 
fective  accelerated  ion  energy  ranges  for  the  4.44  MeV  12C  line 
( white  boxes)  and  the  2.223  MeV  neutron  capture  line  (gray 
boxes)  for  several  power-law  spectral  indexes  and  for  accelerated 
a/proton  =  0.1  (Fig.  21a)  and  0.5  (Fig.  21  b).  For  the  neutron 
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Fig.  22. — Effective  accelerated  ion  energies  for  production  of  neutrons  arriving 
at  Earth  before  ( gray  boxes)  and  after  ( white  boxes)  a  time  4rr  =  850  s  after  the 
release  of  the  accelerated  ions  into  the  magnetic  loop  at  the  Sun  as  a  function  of 
spectral  index.  The  effective  energy  ranges  are  the  50%  yield  range  as  defined  in  the 
text.  Panel  (a)  shows  accelerated  a/proton  =  0.1  and  (b)  shows  a/proton  =  0.5. 
For  these  calculations,  X  —  300,  8  —  0.2,  L  =  1 1,500  km,  and  90 bS  =  75°. 

capture  line,  X  —  300,  8  =  0.2,  3He/H  =  3  x  1 0  5,  and  90 bS  = 
75°.  The  effective  ion  energy  ranges  for  the  4.44  MeV  line  are 
similar  to  those  for  the  6. 129  MeV  160  line  shown  in  Figure  20 
since  the  production  cross  sections  for  the  two  lines  are  similar. 
For  hard  spectra  ( s  <  4),  the  effective  ion  energies  for  the  neu¬ 
tron  capture  line  extend  to  very  high  energies,  up  to  ~  1 00  MeV 
nucleon  ,  and  are  generally  much  higher  than  those  producing 
the  12C  line.  However,  for  soft  spectra  (s  >  4),  the  effective  ion 
energies  for  the  neutron  capture  line  can  be  very  low  (<7  MeV 
nucleon-1)  due  to  neutron  production  by  the  exothermic  ( a,n ) 
reactions  on  heavy  elements.  For  such  soft  spectra  the  effective 
ion  energies  for  the  neutron  capture  line  are  even  less  than  those 
for  the  12C  line.  This  behavior  of  the  separation  of  the  effective  ion 
energies  for  the  two  lines  is  reflected  in  Figure  7,  where  the  slope 
of  the  yield  ratio  curve  is  very  steep  for  hard  spectra  but  then 
flattens  and  even  begins  to  rise  for  very  soft  spectra.  The  very 
broad  effective  ion  energy  distribution  for  the  neutron  capture  line 
at  intermediate  spectral  indexes  reflects  the  nearly  equal  con¬ 
tribution  to  the  line  from  both  low-energy  and  high-energy 
neutron  production  reactions.  We  find  that  there  is  essentially 
no  dependence  of  the  neutron  capture  line  effective  ion  energy 
ranges  on  X  and  8.  As  the  flare  location  moves  from  disk  center 
to  limb,  the  effective  energies  shift  to  slightly  lower  energies 
since  the  neutrons  from  higher  energy  reactions  are  generally  pro¬ 
duced  deeper  and  are  subsequently  more  attenuated  by  Compton 
scattering. 

When  a  spectral  index  is  derived  using  the  ratio  of  these  two 
lines,  the  relevant  ion  energies  cover  a  broad  range  of  energies. 
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Fig.  23. — Effective  accelerated  ion  energies  for  production  of  positrons  as  a 
function  of  spectral  index  for  accelerated  3He/4He  =  0  (white  boxes)  and  1  ( gray 
boxes).  The  effective  energy  ranges  are  the  50%  yield  range  as  defined  in  the  text. 
Panel  (a)  shows  accelerated  a/proton  =  0. 1  and  (b)  shows  a/proton  =  0.5. 
Also  shown  are  the  effective  energies  for  the  4.438  MeV  160  line  (black  boxes) 
from  Fig.  21. 

For  hard  spectra,  the  relevant  ion  energy  range  extends  from  a 
few  MeV  nucleon- 1  up  to  and  greater  than  1 00  MeV  nucleon- 1 . 
For  soft  spectra  (5  >  4),  the  relevant  ion  energy  range  is  much 
narrower,  from  less  than  1  to  a  few  MeV  nucleon- 1 .  For  very  soft 
spectra,  the  lack  of  separation  of  the  effective  ion  energies  pro¬ 
ducing  the  two  lines  reduces  the  sensitivity  of  the  ratio  to  the  spec¬ 
tral  index,  as  seen  in  Figure  7. 

Determining  the  effective  accelerated  ion  energies  for  produc¬ 
tion  of  the  neutrons  escaping  from  the  Sun  and  arriving  at  Earth 
is  more  complicated  than  for  a  7-ray  line  since  there  is  no  simple 
yield.  Rather,  there  is  an  escaping  neutron  spectrum  that  results 
in  a  time  history  of  neutrons  of  varying  energies  arriving  at 
Earth.  We  show  in  Figure  22  the  effective  accelerated  ion  ener¬ 
gies  for  neutrons  arriving  at  Earth  before  ( gray  boxes)  and  after 
(white  boxes)  a  time  tm  =  850  s  after  the  release  of  the  accel¬ 
erated  ions  into  the  magnetic  loop  at  the  Sun.  The  effective  en¬ 
ergies  are  shown  for  several  power-law  spectral  indexes  and 
X  =  300, 8  =  0.2, 3He/H  =  3  x  10-5,  and0Obs  =  75°  and  for  ac¬ 
celerated  a/proton  =  0.1  (Fig.  22 a)  and  0.5  (Fig.  22b).  The  choice 
of  4rr  =  850  s  is  motivated  by  the  analysis  of  the  June  4  flare 
presented  in  §  3.1.  The  figure  shows  that  the  effective  ion  ener¬ 
gies  producing  neutrons  that  arrive  at  the  Earth  before  850  s  are 
very  high,  always  greater  than  100  MeV  nucleon-1.  Neutrons 
arriving  after  850  s  are  produced  by  accelerated  ions  of  lower 
energies. 

In  Figure  23  we  show  the  effective  accelerated  ion  energies  for 
the  production  of  positrons  for  accelerated  3He/4He  =  0  (white 
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TABLE  4 

Measurable  Quantities  and  their  Dependences  on  the  Acceleration  and  Physical  Parameters 


Observable 


Acceleration  Parameter  Physical  Parameter 


Narrow  deexcitation  line  fluences  and  ratios 

Narrow  deexcitation  line  shift  and  shape . 

Narrow  deexcitation  line  time  history . 

Electron  bremsstrahlung  time  history . 

Neutron-capture  line  ftuence . 

Neutron-capture  line  time  history . 

Neutron  ftuence  at  Earth . 

Neutron  arrival  time  history  at  Earth . 

Images . 


alp ,  s 
alp ,  s 

^ion  (0 

<h(?) 
alp,  s 

“ion (/)»  a/p,  s 
a/p,  s 

aio„(0»  alp,  s 


Ambient  abundances 
8,  X,  0obs,  n(h) 

8,  X,  L,  n(h) 

L 

8,  X,  0obs,  n(h),  ambient  3He/H 
8,  X,  dobs,  L,  n(h),  ambient  3He/H 
8,  X,  dobs,  n(h) 

8,  X,  dobs,  L,  n(h) 

8,  dobs 


boxes )  and  1  ( gray  boxes )  and  for  accelerated  ct/proton  =  0.1 
(Fig.  23 a)  and  0.5  (Fig.  23 b).  The  figure  shows  that  for  very  hard 
spectra,  the  effective  ion  energies  can  be  extend  to  more  than  a 
GeV  nucleon-1  because  the  positrons  result  mainly  from  very 
high-energy  reactions  producing  charged  pions.  With  a/proton  = 
0. 1 ,  the  production  is  primarily  from  accelerated  proton  reactions 
and  the  lower  energy  of  the  effective  energy  range  is  several  hun¬ 
dred  MeV  nucleon-1.  When  a/proton  is  0.5,  lower  energy  ac¬ 
celerated  a  particle  and  heavy  element  interactions  producing 
radioactive  nuclei  also  contribute,  extending  the  effective  energy 
range  to  lower  energies.  The  very  low-energy  reactions  of  accel¬ 
erated  3He  do  not  contribute  significantly.  For  softer  spectra,  the 
contribution  from  3He  reactions  can  lower  the  lower  range  of 
the  effective  energies  to  less  than  1  MeV  nucleon-1,  the  lowest 
ion  energies  for  any  of  the  observables  considered  here.  The  up¬ 
per  range  is  from  a  few  MeV  nucleon-1  to  ~15  MeV  nucleon-1 
depending  on  the  3He  abundance. 

Because  the  measured  ratio  of  the  0.511  MeV  positron  anni¬ 
hilation  line  fluence  to  the  fluence  of  a  deexcitation  line,  such  as 
the  4.44  MeV  12C  line,  has  also  been  used  to  derive  spectral  in¬ 
dexes,  for  comparison  we  also  show  in  Figure  23  the  effective  ac¬ 
celerated  ion  energy  ranges  for  the  4.44  MeV  line  from  Figure  2 1 . 
We  see  that  when  spectral  indexes  are  determined  from  this  ratio, 
the  effective  ion  energies  are  from  less  than  an  MeV  nucleon- 1 
to  a  few  MeV  nucleon-1  for  very  soft  spectra  and  from  many 
MeV  nucleon- 1  to  over  a  GeV  nucleon- 1  for  very  hard  spectra. 
We  again  caution,  however,  that  the  effective  energies  shown  in 
this  figure  are  for  the  production  of  positrons.  The  escaping 
0.511  MeV  annihilation  photons  can  be  significantly  attenuated 
if  the  positrons  are  produced  deep  in  the  solar  atmosphere.  Pion 
production  can  be  very  deep  since  the  high-energy  ions  respon¬ 
sible  have  very  long  ranges.  As  a  result,  even  for  very  hard  spec¬ 
tra,  the  annihilation  photons  that  do  escape  can  be  mainly  from 
decay  of  radioactive  nuclei  rather  than  pions.  Such  nuclei  are  pro¬ 
duced  by  ions  with  energies  that  are  much  lower  than  those  indi¬ 
cated  in  the  figure.  High-energy  (>10  MeV)  continuum  emission 
from  neutral  and  charged  pion  decay  is  only  produced  by  accel¬ 
erated  ions  (protons  and  a-particles)  with  energies  greater  than 
~200  MeV  nucleon-1. 

3.  ANALYSIS  PROCEDURE  FOR  HIGH-ENERGY 
FLARE  DATA 

In  the  previous  sections  we  have  explored  how  the  measurable 
quantities  of  high-energy  flare  data  depend  on  the  acceleration 
and  physical  parameters  of  the  transport  and  interaction  model. 
The  measurable  quantities  that  we  will  consider  in  the  flare  anal¬ 
ysis  presented  here  are  deexcitation  line  time-dependent  fluxes, 
time-integrated  fluences  and  fluence  ratios;  deexcitation  line  shapes 
and  centroid  shifts;  time-dependent  electron  bremsstrahlung  flux; 


neutron-capture  line  time-dependent  flux  and  fluence;  escaping 
neutron  time-dependent  flux  at  earth;  and  images  of  the  active 
region  at  various  wavelengths.  We  list  these  quantities  in  Table  4 
along  with  their  most  important  dependences  on  the  acceleration 
and  physical  parameters  (see  Table  1).  Most  of  these  parameters 
have  been  defined  earlier.  The  heliocentric  angle  of  the  flare,  0obs, 
is  provided  by  NOAA/SEC  determinations  based  on  optical  im¬ 
ages.  For  flares  located  near  the  solar  limb,  these  determinations 
can  be  uncertain.  In  addition,  the  flare  optical  location  and  the  site 
of  the  nuclear  reactions  may  not  be  the  same. 

We  see  that  each  of  the  measurable  quantities  can  depend  on 
several  parameters.  Restricting  an  analysis  to  one  or  a  small  sub¬ 
set  of  the  measurable  quantities  can  result  in  parameter  estimates 
that  may  be  inconsistent  with  estimates  derived  from  analysis 
of  other  measurable  quantities.  And,  if  the  effect  of  all  of  the  rel¬ 
evant  parameters  is  properly  taken  into  account,  the  uncertainties 
will  be  necessarily  large;  if  only  a  subset  of  the  parameters  is  con¬ 
sidered,  the  uncertainties  will  be  underestimated.  A  useful  and 
reliable  analysis  must  utilize  as  many  measurable  quantities  as 
possible  to  determine  as  many  parameters  as  possible.  But  be¬ 
cause  of  the  many  interdependences  on  the  parameters,  this  can 
be  quite  challenging  and  making  progress  in  constraining  the  pa¬ 
rameters  may  require  additional  assumptions.  In  the  following  we 
analyze  the  high-energy  flare  data  available  for  the  1991  June  4 
flare  within  the  framework  of  the  loop  transport  model.  We  will 
attempt  to  establish  well-constrained  values  for  all  of  the  param¬ 
eters  and  show  that  the  loop  model  can  account  for  the  wide 
variety  of  observables  available  for  this  flare  with  reasonable 
values  for  them.  The  guiding  principle  of  the  analysis  will  be  to 
first  analyze  those  measurable  quantities  that  depend  only  on  one 
or  two  parameters  and  then,  using  those  parameter  estimates,  pro¬ 
ceed  to  other  observables  that  have  more  complicated  depen¬ 
dences.  The  procedure  is  designed  to  take  best  advantage  of  the 
parameter  dependences  and  may  be  useful  in  analysis  of  other 
flares  with  a  similar  range  of  measurable  quantities. 

3.1.  Analysis  of  the  CGRO  OSSE  1991  June  4  Flare  Data 

In  June  of  1991,  active  region  6659  rotated  onto  the  solar  disk 
and  produced  some  of  the  largest  flares  ever  recorded.  It  first 
appeared  at  the  east  limb  on  June  1  and  produced  an  X12+  flare. 
The  peak  of  this  flare  was  missed  by  the  CGRO  instruments 
because  the  satellite  was  in  the  South  Atlantic  Anomaly  (SAA). 
As  there  was  a  high  probability  for  subsequent  intense  flares,  a 
solar  Target  of  Opportunity  was  declared  for  the  Oriented  Scin¬ 
tillation  Spectrometer  Experiment  (OSSE)  on  CGRO,  which  is 
sensitive  to  7-rays  from  0.04  to  >1 50  MeV  and  neutrons  >10  MeV. 
On  June  4,  AR  6659  produced  a  second  X12+  flare  while  OSSE 
was  viewing  the  Sun.  The  flare  was  one  of  the  most  intense  nu¬ 
clear  7-ray  line  flares  observed  to  date;  most  instruments,  from  the 
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soft  X-ray  to  the  7-ray  bands,  saturated  at  its  peak.  Energetic  par¬ 
ticles  from  the  flare  were  observed  in  interplanetary  space  by  the 
EPAC  experiment  on  Ulysses  (McDonald  et  al.  1995),  GOES  and 
IMPS  (Tylka  et  al.  1996).  The  NOAA  Solar-Geophysical  Data 
Comprehensive  Report  gives  the  Hq  flare  location  as  N30E70. 
This  location  corresponds  to  a  heliocentric  angle  of  ~74.5°  at  the 
time  of  the  flare.  The  flare  was  also  observed  at  1 7, 35  and  80  GHz 
by  radiometers  at  the  Nobeyama  Radio  Observatory  in  Japan 
(Ramaty  et  al.  1994). 

During  the  OSSE  observations,  the  Sun  was  kept  within  0.8° 
of  the  axes  of  two  of  the  detectors,  while  the  remaining  two  de¬ 
tectors  “chopped”;  i.e.,  alternately  pointed  at  and  4.5°  off  the  Sun 
at  131.1  s  intervals.  Full  spectra  were  accumulated  from  each  de¬ 
tector  at  8  s  intervals.  Excellent  measurements  of  both  7-rays  and 
neutrons  were  obtained  during  the  rise,  peak,  and  decay  of  the 
event.  The  decay  was  interrupted  by  satellite  night  but  observa¬ 
tions  were  resumed  at  daylight  of  the  next  orbit  and  additional 
observations  were  obtained.  Nuclear  deexcitation  line  and  neutron- 
capture  line  emission  were  detected  ( Murphy  et  al.  1 997)  through¬ 
out  the  first  orbit  (>1000  s)  and  continued  into  the  second. 
Continuum  7-ray  emission  at  >16  MeV  was  detected  during  the 
peak  of  the  flare.  The  high-energy  emission  decayed  much  more 
rapidly  than  the  line  emission,  falling  to  below  detectability  within 
250  s  after  the  7-ray  peak.  The  spectrum  of  the  high-energy  emis¬ 
sion  showed  no  evidence  for  a  pion-decay  feature  ( Del  Signore 
1995),  implying  that  this  emission  was  predominantly  electron 
bremsstrahlung.  A  strong  neutron  signal  was  detected  following 
the  high-energy  7-ray  emission  and  continued  until  the  Sun  was 
occulted  by  the  Earth.  The  June  4  flare  neutron  data  were  dis¬ 
cussed  in  the  analysis  (Murphy  et  al.  1999)  of  the  1991  June  1 
flare  in  which  neutron  emission  was  also  detected  by  OSSE. 

The  June  4  flare  was  so  intense  that  it  produced  severe  satura¬ 
tion  effects  in  the  three  detectors  directly  viewing  the  Sun.  Data 
from  the  single  off-pointed  detector  were  largely  unaffected  be¬ 
cause  the  hard  X-ray  photons  were  absorbed  in  the  shield  and 
collimator.  During  the  peak,  only  data  from  this  detector  are  used, 
although  the  statistical  significance  is  reduced  and  the  absorption 
and  scattering  prevents  measurement  of  solar  photons  less  than 
~1  MeV. 

We  begin  the  analysis  with  the  parameter  that  is  most  easily 
established,  the  size  of  the  flare  loop,  which  can  be  simply  deter¬ 
mined  from  images  of  the  flare.  Ideally,  images  of  7-ray  emission 
loop  footpoints,  such  as  obtained  by  RHESSI  from  the  2002  Oc¬ 
tober  28  flare,  would  be  used.  Such  images  are  not  available  for 
the  June  4  flare,  but  some  constraints  can  be  obtained  from  im¬ 
ages  at  other  wavelengths.  For  example,  white-light  images  of 
flares  can  be  used  to  identify  potential  loop  footpoints.  Sakurai 
et  al.  (1992)  provide  white-light  images  of  active  region  6659 
during  the  June  4  flare.  The  third  panel  of  their  Figure  2b  shows 
the  white-light  patches  at  3: 40 : 34  UT,  very  near  the  peak  of  the 
June  4  flare  7-ray  line  emission  at  3: 41:07  UT.  If  these  patches 
are  taken  to  be  footpoints  where  ions  are  interacting,  then  the 
possible  corresponding  loop  lengths  range  from  11,500  km  to 
65,000  km  (assuming  a  semicircular  loop)  depending  on  which 
pair  is  assumed  to  be  associated  with  the  loop  producing  the 
7-rays.  In  the  analysis  following  we  will  consider  three  possible 
loop  lengths:  11,500,  34,500  and  65,000  km. 

The  most  direct  measure  of  the  accelerated  a/proton  ratio  is 
the  measured  fluence  ratio  of  the  a-a  line  complex  and  a  narrow 
line  such  as  the  4.438  MeV  line  of  12C  or,  to  improve  statistical 
significance,  the  sum  of  high-FlP  lines.  However,  the  fluence  ra¬ 
tio  also  depends  on  the  accelerated  ion  spectral  index.  The  spec¬ 
tral  index  can  be  determined  from  the  measured  fluence  ratio  of 
the  6.129  MeV  160  and  1.634  MeV  20Ne  lines,  but  this  fluence 


Fig.  24. — Two-parameter  99%  confidence  contours  for  combinations  of  ac¬ 
celerated  a/proton  ratio  and  spectral  index  from  the  a-a  complex-high  FIP  line 
fluence  ratio  comparison  (dotted  curve)  and  the  6.129  MeV  160-1.634  MeV 
20Ne  line  fluence  ratio  comparison  (solid  curve)  for  the  June  4  flare.  The  gray 
area  represents  the  overlap  of  the  99%  confidence  regions  from  both  analyses. 

ratio  also  depends  on  the  accelerated  a/proton  ratio.  Here  we 
will  simultaneously  determine  the  cr/proton  ratio  and  the  spectral 
index  by  calculating  values  of  the  two  fluence  ratios  for  a  range 
of  assumed  cc/proton  ratios  and  spectral  indexes  and  compar¬ 
ing  them  with  the  measured  fluence  ratios.  The  measured  June  4 
fluence  ratio  of  the  a-a  line  complex  and  the  sum  of  high-FIP 
narrow  lines  obtained  with  Sun-pointed  detectors  late  in  the  flare 
was  0.66  ±  0.04  (Murphy  et  al.  1997).  The  flare-averaged  mea¬ 
sured  fluence  ratio  of  the  6.129  MeV  160  and  1.634  MeV 
20Ne  lines  integrated  over  the  entire  flare  0.66  ±  0.07  (see  er¬ 
ratum  to  Murphy  et  al.  1997).  We  show  in  Figure  24  the  99% 
confidence  contours  from  each  of  these  comparisons,  with  the 
gray  area  representing  the  overlapping  99%  confidence  region. 
We  see  that  the  accelerated  a/proton  ratio  was  at  least  0.4  and 
the  spectral  index  is  3.9  ±  0.2.  Murphy  et  al.  (1997)  also  found 
that  the  a/proton  ratio  was  closer  to  0.5  than  0.1.  While  con¬ 
sistent  within  the  uncertainties  with  the  index  (4.4  ±  0.3)  deter¬ 
mined  by  Murphy  et  al.  (1997),  the  index  derived  here  is  harder 
due  to  the  updated  cross  sections  used  here.  We  see  from  Fig¬ 
ure  20  that  the  ion  energies  relevant  to  a  spectral  index  determined 
from  the  O/Ne  line  ratio  is  from  ~1  to  several  MeV  nucleon-1. 
Murphy  et  al.  (1997)  showed  that  within  the  measurement  un¬ 
certainties,  there  was  no  evidence  for  variation  of  the  spectral  in¬ 
dex  with  time.  We  therefore  assume  a  constant  power- law  spectral 
index. 

The  parameters  determined  thus  far  do  not  depend  on  the  at¬ 
mospheric  model  n(h).  Because  the  atmospheric  model  cannot  be 
characterized  by  a  simple  parameterization,  we  will  carry  out  the 
analysis  procedure  first  assuming  the  Avrett  atmospheric  model 
representing  a  sunspot  and  then  repeat  the  procedure  assuming  the 
RHESSI/ Avrett  model  representing  possible  flare  conditions  and 
compare  the  results. 

3.1.1.  Analysis  Assuming  the  Avrett  Sunspot  Atmospheric  Model 

The  measurable  quantity  that  depends  on  the  fewest  yet- 
undetermined  parameters  is  the  deexcitation  line  centroid  shift. 
As  discussed  in  §  2.4.2,  this  shift  depends  primarily  on  the  inter¬ 
acting  angular  distribution,  which  in  turn  depends  on  the  mag¬ 
netic  convergence  (6)  and  level  of  PAS  (1).  The  energy  of  the 
photon  emitted  by  deexcitation  from  the  first  excited  state  of  12C 
in  its  rest  frame  is  4.438  MeV.  The  line  centroid  during  the  flare 
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Fig.  25. — Two-parameter  confidence  contours  for  combinations  of  <5  and  X 
resulting  from  the  line  shift  analysis  of  the  June  4  flare. 

was  measured  to  be  4.442  ±  0.005  MeV  by  fitting  a  Gaussian 
profile  to  the  spectral  data.  We  compare  this  measured  centroid 
with  predicted  centroids  calculated  for  angular  distributions  aris¬ 
ing  from  various  combinations  of  2  and  6,  accelerated  a/proton  = 
0.5,  and  the  range  of  spectral  index  determined  above.  The  pre¬ 
dicted  centroids  were  obtained  by  fitting  a  Gaussian  profile  to 
calculated  line  profiles  (e.g.,  see  Fig.  9).  For  each  combination  of 
6  and  2  (ands),  we  calculate  A  <r;  i.e.,  how  many  1  a  uncertainties 
the  measured  line  centroid  differs  from  the  calculated  centroid.  In 
Figure  25,  we  show  contours  of  A  a  =  2.3,  4.6,  and  9.2  corre¬ 
sponding  to  68%,  90%,  and  99%  confidence  (for  2  parameters  of 
interest)  that  the  calculated  shift  is  not  consistent  with  the  mea¬ 
sured  shift.  While  simultaneously  low  values  of  6  and  X  can  be 
rejected  (corresponding  to  downward-directed  angular  distribu¬ 
tions),  there  is  still  a  wide  range  of  <5-2  combinations  that  is  al¬ 
lowed.  Because  of  the  location  of  this  flare  (74.5°  heliocentric 
angle),  the  shift  of  the  line  centroid  from  a  loop  perpendicular  to 
the  solar  surface,  as  assumed  here,  is  not  sensitive  to  the  interac¬ 
ting  particle  angular  distribution  (see  §  2.4.2)  and  so  does  not 
provide  strong  constraints  on  6  and  X. 

A  different  and  more  restrictive  constraint  on  6  and  2  can  be 
obtained  by  considering  the  time  history  of  the  prompt  deexci¬ 
tation  lines.  Within  our  loop  transport  model,  a  necessary  com¬ 
ponent  of  this  analysis  is  the  time  history  of  the  release  of  the 
accelerated  ions  at  the  top  of  loop,  al0n(t),  which  is  not  known. 
We  will  assume  two  cases:  (1)  the  case  in  which  the  ion  release  is 
instantaneous  and  the  time  history  of  the  prompt  line  emission  is 
due  entirely  to  transport  and  (2)  the  case  in  which  the  ion  release 
is  extended  in  time.  For  the  instantaneous  release  case  we  will 
assume  two  possible  instantaneous  release  times:  (1)  at  the  be¬ 
ginning  of  the  >300  keV  electron  bremsstrahlung  as  measured 
by  BATSE  (3:40:50  UT)  and  (2)  at  the  beginning  of  the  prompt 
deexcitation  line  emission  (3:41:00  UT).  For  the  extended  re¬ 
lease  case,  we  will  assume  that  aion(0  is  identical  to  the  relativ¬ 
istic  electron  acceleration  release  time,  ae(t),  and  further  assume 
that  ae(t)  is  identical  to  the  time  history  of  the  relativistic  electron 
bremsstrahlung  flux,  UbrcmU)-  Because  relativistic  electron  trans¬ 
port  should  be  very  fast,  the  bremsstrahlung  emission  time  history 
should  be  very  nearly  identical  to  ae(t).  For  each  release  scenario 
we  calculate  predicted  prompt  deexcitation  line  time  histories  for 
various  combinations  of  <5  and  2  and  compare  them  with  the  time 
history  of  the  OSSE  measured  summed  high-FIP  narrow  y-ray 
lines  for  each  of  the  three  loop  lengths. 


Fig.  26. — Comparison  of  measured  high-FIP  deexcitation  line  time  history 
at  the  peak  of  the  June  4  flare  with  calculated  time  histories  for  various  values 
of  PAS. 


For  the  instantaneous  release-time  scenarios,  we  found  that 
the  best  fit  was  achieved  for  a  loop  length  of  34,500  km.  Loop 
lengths  of  1 1 ,500  and  65,000  km  did  not  provide  acceptable  fits 
for  either  of  the  release  times.  For  the  34,500  km  loop  length, 
only  the  release  time  at  3 : 41 : 00  UT  provided  an  acceptable  fit. 
The  3:40:50  UT  release  time  was  too  early  so  that  even  with  6-2 
combinations  providing  maximum  delay  of  the  emission,  trans¬ 
port  could  not  provide  an  adequate  fit  to  the  time  history.  We 
show  in  Figure  26  the  OSSE  summed  prompt  high-FIP  deexcita¬ 
tion  line  flux  time  history  (data  points)  along  with  time  histories 
calculated  for  6  =  0.20  and  several  values  for  PAS:  2  =  20  (sat¬ 
urated),  300, 700, 5000  and  none  (2  — >  oo).  We  show  in  Figure  27 
the  68%,  90%,  and  99%  6-2  confidence  contours  ( thick  curves ) 
obtained  from  this  fit.  The  figure  shows  that  as  the  level  of  mag¬ 
netic  convergence  increases  (resulting  in  a  more  extended  time 
history),  the  level  of  PAS  must  increase  (shortening  the  time  his¬ 
tory).  Regions  of  low  convergence-high  PAS  and  high  conver¬ 
gence-low  PAS  are  not  acceptable.  We  also  show  in  Figure  27 


Fig.  27. — Two-parameter  confidence  contours  ( thick  curves)  for  combina¬ 
tions  of  8  and  X  resulting  from  the  time  history  analysis  of  the  June  4  flare 
assuming  instantaneous  accelerated  ion  release  at  3 : 41 : 00  UT.  Also  shown  are 
the  contours  shown  in  Fig.  19  from  the  line  shift  analysis  (thin  curves).  The  gray 
area  represents  the  overlap  of  the  99%  confidence  regions  from  both  analyses. 
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Fig.  28. — Time  histories  of  the  high-FIP  narrow  deexcitation  line  flux  ob¬ 
served  with  OSSE  from  the  June  4  flare  and  the  >300  keV  bremsstrahlung  flux 
observed  with  BATSE.  The  upper  B  ATSE  data  are  without  background  subtrac¬ 
tion  and  the  lower  data  are  with  a  simple  linear  background  subtraction. 

the  contours  from  the  line  shift  analysis  (thin  curves,  see  Fig.  25), 
with  the  gray  region  representing  the  overlapping  99%  confidence 
region. 

We  next  consider  the  case  using  the  bremsstrahlung  time  his¬ 
tory.  We  show  in  Figure  28  the  time  history  of  the  summed  prompt 
high-FIP  narrow  7-ray  lines,  </>highFip(0  ( histogram ),  and  the  time 
history  of  the  >300  keV  bremsstrahlung  emission,  </>brem(0>  mea¬ 
sured  by  BATSE,  where  the  upper  curve  is  without  background 
subtraction  and  the  lower  curve  was  obtained  with  a  simple,  linear 
background  subtraction.  For  comparison,  the  7-ray  and  brems¬ 
strahlung  emissions  have  been  normalized  to  the  same  flux  at  the 
peak  of  the  emission.  N  ear  the  peak,  the  two  emissions  have  a  sim¬ 
ilar  time  history.  At  later  times,  while  the  two  time  histories  are 
similar,  the  line  emission  is  stronger  than  the  bremsstrahlung,  as 
noted  by  Murphy  et  al.  (1997).  This  may  suggest  that,  while  later 
in  the  flare  their  acceleration  processes  may  differ,  at  the  peak  of 
the  flare  both  the  ions  and  electrons  are  similarly  accelerated.  We 
therefore  assume  that  the  electron  bremsstrahlung  time  history 
also  represents  the  ion  acceleration  release  time  history  during  the 
first  peak  and  we  concentrate  this  part  of  the  analysis  on  this  peak 
only,  from  3:40:50  to  3:42:15  UT. 

For  this  case  we  find  that  all  loop  lengths  considered  could 
provide  adequate  fits.  However,  even  the  99%  confidence  re¬ 
gions  for  the  34,500  and  65,000  km  loop  lengths  did  not  overlap 
with  the  99%  confidence  region  obtained  from  the  line  shift 
analysis.  In  addition,  these  loop  lengths  only  allowed  values  of  6 
less  than  about  0.10.  Such  small  values  of  6  do  not  provide  the 
increase  of  magnetic  field  from  the  corona  to  the  photosphere 
expected  to  exist  in  flare  loops,  so  we  reject  these  longer  loop 
lengths  as  not  physically  acceptable.  Only  the  11,500  km  loop 
length  provides  a  region  of  acceptable  b  and  7-values.  We  show 
in  Figure  29  the  68%,  90%,  and  99%  8-X  confidence  contours 
(thick  curves)  obtained  from  this  fit.  We  also  show  the  contours 
from  the  line  shift  analysis  (thin  curves,  see  Fig.  25),  with  the  gray 
region  representing  the  overlapping  99%  confidence  region. 

Both  of  the  two  best-fitting  scenarios  (instantaneous  produc¬ 
tion  with  L  =  34,  500  km  and  BATSE  time  history  with  L  = 
11,500  km)  provided  adequate  fits  to  the  measured  7-ray  line 
time  history.  Longer  loop  lengths  either  provided  unacceptable 
fits  or  resulted  in  parameter  values  that  were  not  consistent  with 
those  obtained  from  the  line  shift  analysis  or  were  physically  un¬ 
reasonable,  implying  that  the  June  4  loops  were  relatively  com¬ 
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Fig.  29. — Two-parameter  confidence  contours  (thick  curves)  for  combina¬ 
tions  of  8  and  2  resulting  from  the  time  history  analysis  of  the  June  4  flare  as¬ 
suming  the  accelerated  ion  release  time  history  is  given  by  the  >300  keV  BATSE 
time  history.  Also  shown  are  the  contours  shown  in  Fig.  19  from  the  line  shift 
analysis  (thin  curves).  The  gray  area  represents  the  overlap  of  the  99%  confidence 
regions  from  both  analyses. 

pact,  less  than  ^60,000  km.  If  the  ion  release  time  history  were 
instantaneous,  the  loop  length  must  be  of  intermediate  length, 
~34,000  km.  If  the  release  time  history  were  more  extended  sim¬ 
ilar  to  the  BATSE  time  history,  the  loop  must  be  more  compact, 
~1 1,000  km. 

The  99%  confidence  regions  of  the  b- 7  parameter  space  for  the 
two  scenarios  (see  Figs.  27  and  29)  are  similarly  shaped  and  po¬ 
sitioned  in  the  6-7  space  so  that  using  values  of  6  and  7  drawn 
from  either  region  will  produce  similar  calculated  results  for 
measurable  quantities  that  depend  on  these  parameters.  We  will 
next  analyze  the  neutron  capture  line,  and  from  §  2.4.3  we  know 
that  its  time  history  does  not  depend  significantly  on  loop  length 
forL  <  80,000  km  so,  again,  either  scenario  will  suffice.  We  will 
use  the  loop  length  obtained  from  the  BATSE  time  history  since 
we  believe  the  ion  acceleration  time  history  would  more  likely 
be  similar  to  it  rather  than  to  an  instantaneous  injection.  For  the 
remaining  analyses  we  will  assume  a  range  of  8- 7  combinations 
from  the  acceptable  region  obtained  with  the  BATSE  profile  and 
the  11,500  km  loop  length. 

We  next  analyze  the  neutron-capture  line,  considering  both  the 
shape  of  its  time  history  and  its  absolute  fluence  relative  to  the 
other  emissions.  From  Table  4  all  of  the  parameters  that  this  line 
depends  on  have  now  been  constrained  by  our  procedure  except 
the  ambient 3  He/H  ratio.  We  calculate  predicted  neutron  capture¬ 
line  time  histories  by  convolving  the  time  histories  calculated  for 
instantaneous  release  (see  §  2. 4. 3. 2)  with  the  ion  acceleration 
release  time  history,  ajon(t).  Again,  amn(t)  is  not  known,  but  we 
cannot  use  the  electron  bremsstrahlung  time  history  as  we  have 
in  the  previous  analysis  of  the  deexcitation  line  time  history  at 
the  peak  of  emission  because  we  require  the  time  history  through¬ 
out  the  entire  flare,  and  we  know  that  at  later  times  the  electrons 
and  ions  are  not  similarly  accelerated  (see  Fig.  28).  Instead,  we 
will  use  the  measured  high-FIP  deexcitation  line  time  history, 
<AhighFip(0>  to  represent  the  acceleration  release  time  history. 
Deexcitation  lines  are  prompt,  and  their  time  history  is  essentially 
identical  to  the  nuclear  interaction  time  history.  Ion  transport  times 
from  release  to  interaction  are  generally  much  shorter  than  neutron 
transport  times  from  production  to  capture  by  H  (compare  Figs.  2 
and  17).  As  a  result,  the  capture  line  time  history  should  not  be 
sensitive  to  the  small  time  delay  between  accelerated  ion  release 
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Time  (UT  on  4  June  1991) 

Fig.  30. — Comparison  of  the  2.223  MeV  neutron  capture  line  time  history 
measured  with  OSSE  from  the  June  4  flare  with  a  time  history  calculated  with 
the  parameters  indicated  in  the  text.  Also  shown  is  the  time  history  of  the  high- 
FIP  narrow  deexcitation  line  flux  assumed  to  represent  the  accelerated  ion 
release  time  history.  The  dotted  curves  represent  the  ±  1  o  uncertainties  for  the 
calculated  curve. 

and  the  nuclear  interaction  producing  the  neutron.  For  the  neutron- 
capture  line  analysis,  therefore,  using  the  prompt  deexcitation  line 
time  history  to  represent  the  ion  release  time  history  is  adequate. 
The  absolute  flux  of  the  deexcitation  lines  fixes  the  normalization 
of  the  predicted  neutron  capture  line  time  history. 

We  use  the  neutron  production  code  to  calculate  expected 
neutron-capture  line  time  histories  at  the  heliocentric  angle  of  the 
June  4  flare  for  a  range  of  ambient  3He/H.  We  use  a  range  of  spec¬ 
tral  indexes  drawn  from  the  acceptable  range  of  3.9  ±  0.2  estab¬ 
lished  from  the  combined  a-a/high  F1P  line  ratio  and  O/Ne  line 
ratio  analyses  above  and  assume  the  power  law  extends  unbro¬ 
ken  to  10  GeV.  We  use  a  range  of  8-2  combinations  drawn  from 
the  acceptable  confidence  regions  established  from  the  previous 
line  centroid  shift  and  line  time  history  analyses  (see  Figs.  27  and 
29).  We  compare  the  calculated  time  histories  with  the  measured 
(Murphy  et  al.  1997)  time  history,  using  \2  as  the  measure  of 
goodness  of  fit.  An  example  of  such  a  fit  is  shown  in  Figure  30, 
where  the  solid  curve  is  the  calculated  time  history  for  s  =  4,  sat¬ 
urated  PAS,  8  =  0.35  and  3He/H  =  1  x  10-5.  The  dotted  curves 
represent  the  ±  1  a  uncertainties  for  the  calculated  curve.  Unfor¬ 
tunately,  the  allowed  ranges  for  8,  X,  and  s  we  have  established 
previously  do  not  allow  a  well-determined  value  for  3He/H.  We 
find  that  the  fits  to  the  neutron-capture  line  time  history  could 
only  provide  a  3He/H  upper  limit  of  4.5  x  10~5  (3  a). 

To  compare  our  3He/H  upper  limit  with  values  obtained  with 
other  techniques,  we  use  a  determination  of  photospheric  3He/4He 
obtained  from  recent  solar  wind  measurements  (Geiss  &  Gloeckler 
1998).  Taking  into  account  the  effects  of  separation  and  fraction¬ 
ation,  they  derived  the  present  value  to  be  3He/4He  =  (3.8  ± 
0.5)  x  10~4.  This  value  is  consistent  with  the  value  that  can  be 
derived  from  the  measurement  by  Mahafiy  et  al.  (1998)  of  D/H 
and  3He/4He  in  the  Jovian  atmosphere  using  the  Galileo  probe. 
The  Jovian  abundances  are  approximately  those  of  the  protosolar 
cloud.  In  the  Sun,  all  of  the  protosolar  deuterium  has  been  con¬ 
verted  to  3He.  Taking  into  account  the  small  additional  produc¬ 
tion  of  3He  over  the  solar  history  and  the  small  decrease  of  4He/ FI 
due  to  settling  of  4He  (Geiss  &  Gloeckler  1998),  we  find  that 
this  estimate  of  the  present  value  of  3He/4He  is  consistent  with 
the  above  value  from  the  solar  wind  measurements.  To  convert  to 
3He/H,  the  value  of  3He/4He  must  be  multiplied  by  the  4He/H 
ratio  in  the  photosphere  of  the  ambient,  flaring  plasma.  We  use  the 
value  of  0.097  from  Grevesse  et  al.  (1996).  Because  we  are  dealing 


with  flare  plasma  where  variations  of  4He/H  have  been  observed, 
we  adopt  an  uncertainty  of  ±0.02  for  this  ratio  that  reflects  the 
range  of  values  obtained  by  a  number  of  investigators  (e.g.,  see 
summary  by  Feldman  et  al.  2005).  This  gives  3He/H  =  (3.7  ± 
0.9)  x  10~5  for  our  best  estimate  of  the  photospheric  value.  Our 
upper  limit  of  3Fle/H  <  4.5  x  10~5  derived  from  the  7-ray  mea¬ 
surements  is  consistent  with  this  value. 

We  can  now  use  the  current  best-estimate  3He/H  value  of  (3 .7  ± 
0.9)  x  10  3  to  provide  a  better  constraint  on  the  spectral  index  s 
for  this  flare.  Using  this  range  of  3He/H  values  we  vary  the  spec¬ 
tral  index  and  use  8-2.  combinations  drawn  from  the  acceptable 
ranges  established  above  to  provide  the  best  fit  to  the  neutron- 
capture  line  time  history.  We  find  that  good  fits  are  achieved  with 
3.93  <  s  <  4.05  (±1  a).  This  allowable  range  lies  within  the 
allowable  range  of3.7<s<4.1  obtained  from  the  combined 
a-a/high  FIP  ratio  and  O/Ne  ratio  analyses  and  provides  our  best 
estimate  for  the  index:  3.93  <  s  <  4.05  or  3.99  ±  0.06  (±1  a). 

Because  no  evidence  of  a  pion-decay  emission  feature  was 
seen  in  the  OSSE  >16  MeV  7-ray  spectrum  (Del  Signore  1995), 
the  high-energy  emission  (  detectable  only  during  the  peak  of  the 
7-ray  line  emission)  was  attributed  to  relativistic  electron  brems- 
strahlung.  Using  the  calculations  of  Murphy  &  Ramaty  (1984) 
and  Murphy  et  al.  (1987)  and  assuming  that  the  accelerated-ion 
power-law  spectrum  continues  unbroken  to  energies  >100  MeV 
with  the  spectral  index  range  determined  above,  we  estimate  the 
total  7-ray  fluence  at  100  MeV  due  to  pion-decay  during  the  first 
interval  would  be  ~0.03  photons  MeV-1  cm~2.  The  measured 
(Del  Signore  1995)  fluence  at  100  MeV  was  0.12  photons  MeV-1 
cnr2.  We  therefore  confirm  that  the  bulk  of  the  high-energy  emis¬ 
sion  was  electron  bremsstrahlung.  If  the  accelerated-ion  spectrum 
were  to  steepen  above  100  MeV,  the  expected  pion-decay  fluence 
would  be  less  (see  the  neutron  analysis  below). 

Since  all  of  the  parameters  listed  in  Table  4  have  now  been 
determined,  the  remaining  observables  (the  neutron  arrival  time 
history  and  absolute  fluence)  will  provide  a  consistency  check  on 
those  parameters.  In  addition,  since  the  neutrons  have  a  different 
parameter  sensitivity  than  the  other  measurable  quantities,  their 
measurement  can  instruct  how  to  modify  the  other  parameter  es¬ 
timates  and  iterate  the  above  procedure  to  provide  better  param¬ 
eter  constraints.  We  use  the  neutron  production  code  to  calculate 
escaping  neutron  kinetic-energy  spectra  for  the  heliocentric  an¬ 
gle  of  the  June  4  flare,  0obs  =  74.5°.  For  the  ion  acceleration 
release  time  history,  amn{t),  we  again  use  the  prompt  deexcitation 
line  time  history,  0highFip(O-  We  propagate  the  escaping  neu¬ 
trons  to  Earth,  taking  into  account  decay  of  the  neutrons,  and 
convolve  the  arriving  neutron  spectra  with  the  OSSE  neutron 
response  (see  Appendix),  producing  predicted  neutron  energy- 
loss  count  spectra.  These  spectra  are  integrated  above  23  MeV 
(OSSE  channel  2)  and  compared  with  the  corresponding  mea¬ 
sured  count  rates.  The  statistical  counting  rate  uncertainties  have 
been  increased  by  adding  in  quadrature  an  additional  uncertainty 
of  30%  to  account  for  uncertainties  in  the  OSSE  neutron  response 
function  (see  the  Appendix). 

For  all  of  the  8-2  combinations  allowed  from  the  previous 
analyses  (see  Figs.  27  and  29)  and  the  index  range  of  3 .99  ±  0.06 
determined  above,  we  find  that  the  number  of  neutrons  arriving  at 
Earth  is  overpredicted,  with  the  discrepancy  being  larger  at  early 
times.  This  can  be  seen  in  Figure  3 1  comparing  the  measured  neu¬ 
tron  count  rate  with  a  predicted  rate  {solid  curve )  calculated  for 
saturated  PAS,  8  =  0.35  and  s  at  its  upper  (i.e.,  steeper  spectrum) 
1  a  limit  of  4.05.  The  neutron  count  rate  is  overpredicted  at  times 
less  than  about  850  s  after  the  peak  of  the  deexcitation  line 
emission.  Figure  22  shows  that  for  s  =  4,  neutrons  arriving 
before  850  s  are  produced  by  accelerated  ions  with  energies 
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Fig.  3 1 . — Comparison  of  the  >23  MeV  neutron  count  rate  time  history  mea¬ 
sured  with  OSSE  from  the  June  4  flare  with  count  rate  time  histories  calculated 
with  the  parameters  indicated  in  the  text.  The  solid  curve  is  for  an  unbroken 
power-law  accelerated-particle  spectrum  and  the  dotted  curve  is  for  a  spectrum 
cutoff  at  Ec  =  125  MeV  nucleon-1.  Also  shown  is  the  time  history  of  the  high- 
FIP  narrow  deexcitation  line  flux  {histogram)  assumed  to  represent  the  accel¬ 
erated  ion  release  time  history. 

>200  MeV  nucleon- 1 .  Reducing  their  number  relative  to  lower 
energy  neutrons  can  be  achieved  if  the  accelerated-particle  spec¬ 
trum  steepens  with  energy,  as  long  as  this  does  not  affect  the  other 
measurable  quantities  previously  analyzed.  This  is  possible  since 
deexcitation  lines  and  the  neutrons  responsible  for  the  neutron 
capture  line  are  produced  by  accelerated  ions  with  lower  energies 
than  the  ions  producing  the  neutrons  escaping  into  space  and  sur¬ 
viving  to  Earth  (see  §  2.4.4  ).  Such  steepening  has  been  observed  in 
SEP  spectra  from  both  impulsive  and  gradual  flares  (Tylka  et  al. 
2005, 2002;  Mason  et  al.  2002).  Here  we  will  simulate  such  steep¬ 
ening  by  introducing  a  sharp  cutoff  of  the  power-law  spectrum  at 
an  energy  Ec,  realizing  that  an  actual  particle  spectrum  most  likely 
would  not  cut  off  suddenly  but  would  steepen  gradually  with  en¬ 
ergy.  We  also  show  in  Figure  3 1  ( dotted  curve )  the  predicted  neu¬ 
tron  arrival  time  history  calculated  for  the  same  spectral  index  s  = 
4.05  but  with  a  cutoff  energy  Ec  =  125  MeV  nucleon-1.  We  have 
confirmed  that  this  spectrum  still  predicts  a  neutron-capture  line 
time  history  that  provides  an  acceptable  fit  to  the  data  and  that  the 
other  calculated  measurable  quantities  are  not  affected.  As  dis¬ 
cussed  in  §  2.4.4,  the  other  observables  are  primarily  produced  by 
accelerated  ions  with  energies  less  than  125  MeV  nucleon-1. 

3.1.2.  Analysis  Assuming  the  RHESSI/Avrett  Atmospheric  Model 

Starting  with  the  previously  determined  parameters  that  are 
independent  of  the  atmospheric  model  (s  =  3.9  ±  0.2,  and 
a/proton  >  0.4),  we  repeat  the  previous  analysis  but  using  the 
RHESSI/Avrett  atmospheric  model.  As  for  the  Avrett  analysis,  we 
begin  by  attempting  to  constrain  8  and  A.  As  discussed  in  §  2.3.2, 
the  interacting  particle  angular  distribution  (and  thus  the  shift  of 
the  deexcitation  line  centroid)  is  similar  for  both  models  of  the 
atmosphere,  except  for  very  low  levels  of  PAS,  where  somewhat 
greater  redshifts  result  from  the  slightly  more  downward-directed 
interacting  angular  distribution.  As  for  the  Avrett  atmosphere,  the 
line-shift  analysis  does  not  provide  any  significant  constraint  on  <5 
and  A  for  this  flare  location. 

We  now  attempt  to  constrain  S  and  A  by  considering  the  prompt 
deexcitation  line  time  history  using  the  same  ion  acceleration  re¬ 
lease  time  histories  assumed  in  the  Avrett  analysis.  Because  the 
ion  interaction  time  histories  for  this  atmosphere  are  much  shorter 
than  for  the  Avrett  atmosphere  (see  §  2.3 . 1),  we  found  that  the  time 
histories  for  neither  of  the  instantaneous  release  times  could  pro¬ 


vide  acceptable  fits;  no  combination  of  8  and  A  could  produce 
a  time  history  extended  enough  to  fit  the  line  time  history.  If  the 
RHESSI/Avrett  atmosphere  is  applicable  during  this  flare,  the  im¬ 
plication  of  this  analysis  is  that  accelerated-particle  release  cannot 
be  instantaneous.  Using  the  electron  bremsstrahlung  time  history 
as  the  ion  acceleration  release  time  history,  we  again  find  that  all 
of  the  loop  lengths  could  provide  adequate  fits.  But,  as  discussed 
in  §2.3.1  and  shown  in  Figure  2 \f  the  dependence  of  the  time  his¬ 
tory  on  6  and  A  is  weaker  for  the  RHESSI/Avrett  atmosphere  than 
for  the  Avrett  atmosphere.  This,  combined  with  the  relatively  poor 
statistical  quality  of  the  OSSE  off-pointed  data,  resulted  in  no  sig¬ 
nificant  constraint  on  6  and  A. 

With  no  constraints  on  8  and  A  and  the  ±0.2  uncertainty  on 
the  spectral  index,  no  significant  constraint  on  the  photospheric 
3He/H  ratio  can  be  established  using  the  time  history  of  the  neu¬ 
tron  capture  line.  We  now  fix  the  3He/H  ratio  at  the  independently 
determined  value  of  (3.7  ±  0.9)  x  10-5  and  reverse  the  analysis 
attempting  to  place  constraints  on  8  and  A.  But,  except  for  exclu¬ 
sion  of  8-A  combinations,  where  both  parameters  are  very  small 
(corresponding  to  strongly  downward-directed  interacting  ion 
angular  distributions),  only  very  broad  allowable  ranges  of  8- A 
combinations  can  be  established. 

For  the  neutron  arrival  time  history,  we  find  that  a  good  fit  of 
the  predicted  count  rate  to  the  data  could  not  be  obtained,  regard¬ 
less  of  the  A -A  combination  assumed.  Similar  to  the  result  for  the 
Avrett  atmosphere,  the  count  rate  was  overpredicted  at  early  times. 
We  again  find  that  introducing  a  cutoff  of  the  accelerated  ion 
spectrum  above  ^125  MeV  nucleon-1  provides  a  good  fit  with¬ 
out  affecting  the  other  measurable  quantities. 

4.  DISCUSSION 

Using  recently  updated  cross  sections  for  deexcitation  7-ray 
lines  (Kozlovsky  et  al.  2002)  and  neutron  production  (Hua  et  al. 
2002),  we  have  calculated  a  number  of  measurable  quantities 
associated  with  narrow  7-ray  deexcitation  lines,  the  neutron  cap¬ 
ture  line,  and  escaping  neutron  spectra.  All  of  these  calculations 
have  been  performed  within  the  context  of  a  magnetic  loop  model 
that  provides  a  physical  basis  for  the  interacting  accelerated  ion 
time  delays  and  anisotropic  angular  distributions  that  previously 
have  only  been  assumed.  The  model  includes  energy  losses  due  to 
Coulomb  collisions,  removal  by  nuclear  reactions,  magnetic  mir¬ 
roring  in  the  convergent  flux  tube,  and  MHD  pitch-angle  scatter¬ 
ing  ( PAS)  in  the  corona.  The  measurable  quantities  depend  on  the 
parameters  of  the  model  (see  Table  1),  which  are  of  two  types: 
acceleration  and  physical.  Acceleration  parameters  are  associated 
with  the  acceleration  processes:  accelerated  particle  power-law 
spectral  index  (5),  acceleration  release  time  history  [aion(f)],  and 
accelerated  particle  elemental  composition.  Physical  parameters 
are  associated  with  the  conditions  of  the  magnetic  loop  within 
which  the  accelerated  particles  transport  and  interact  with  the 
ambient  atmosphere:  level  of  PAS  (A),  magnetic  field  conver¬ 
gence  (8),  loop  length  ( L ),  and  the  atmospheric  elemental  com¬ 
position  and  density  and  temperature  height  distributions  [«(/?) 
and  T(h)]. 

The  measurable  quantities  we  considered  were:  narrow  de¬ 
excitation  line  shapes,  shifts,  fluences,  and  time  histories;  neutron 
capture  line  fluence  and  time  history;  and  escaping  neutron  Ali¬ 
enee  and  time  history.  Each  of  these  measurable  quantities  de¬ 
pends  to  varying  degree  on  some  subset  of  the  acceleration  and 
physical  parameters  (see  Table  4).  We  explored  these  dependences 
in  detail,  presenting  a  number  of  figures  expressing  the  most  sig¬ 
nificant  effects  for  each  observable.  To  present  these  dependences 
for  every  possible  combination  of  physical  and  acceleration  pa¬ 
rameters  would  be  impossible.  Instead,  we  carefully  chose  ranges 
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of  parameter  values  that  most  effectively  convey  these  often 
complex  interdependences. 

We  find  that  deexcitation  line  yields  are  not  sensitive  to  phys¬ 
ical  parameters  other  than  the  ambient  composition,  being  most 
sensitive  to  the  acceleration  parameters.  The  line  shapes  and 
shifts,  however,  are,  since  the  physical  parameters  affect  the  an¬ 
gular  distribution  of  the  interacting  accelerated  particles.  Phys¬ 
ical  parameters  also  affect  the  delay  between  accelerated  ion  re¬ 
lease  and  interaction,  which  affects  the  line  time  history.  Similarly, 
total  neutron  production  is  also  more  sensitive  to  the  acceleration 
parameters.  However,  the  number  and  spectrum  of  neutrons  escap¬ 
ing  from  the  Sun  along  a  given  direction  and  the  neutron  capture 
line  fluence  and  time  history  are  very  sensitive  to  the  physical  param¬ 
eters,  again  due  to  the  effect  of  the  interacting  ion  angular  distribu¬ 
tion.  Increased  production  of  neutrons  in  interactions  of  accelerated 
ions  moving  downward  toward  the  Sun  increases  the  neutron- 
capture  line  yield  but  decreases  the  yield  of  escaping  neutrons. 

Because  each  of  the  measurable  quantities  can  depend  on  some 
or  all  of  the  same  acceleration  or  physical  parameters,  deriving 
well-constrained  values  for  those  parameters  can  be  difficult.  Re¬ 
stricting  an  analysis  to  one  or  a  small  subset  of  the  measurable 
quantities  can  result  in  parameter  estimates  that  may  be  incon¬ 
sistent  with  estimates  derived  from  analysis  of  other  measurable 
quantities  and  incorrectly  estimated  uncertainties.  Acknowledg¬ 
ing  that  a  meaningful  analysis  must  utilize  as  many  measurable 
quantities  as  possible,  we  applied  the  calculations  in  an  analysis  of 
the  7-ray  and  neutron  observations  of  the  1991  June  4  solar  flare 
obtained  with  CGRO  OSSE,  one  of  the  best-observed  high-energy 
flares.  The  analysis  technique  applied  to  this  flare  provides  an 
algorithm  that  can  be  applied  to  other  flares  with  measurements 
covering  a  similar  wide  range  of  measurable  quantities.  The  guid¬ 
ing  principle  of  the  analysis  is  to  first  analyze  those  measurable 
quantities  that  depend  only  on  one  or  two  parameters  and  then, 
using  those  parameter  estimates,  proceed  to  other  observables  that 
have  more  complicated  dependences.  The  procedure  is  designed 
to  take  best  advantage  of  the  parameter  dependences.  We  showed 
that  the  loop  model  can  account  for  the  wide  variety  of  ob¬ 
servables  available  for  this  flare  with  reasonable  values  for  the 
parameters. 

We  began  the  analysis  by  obtaining  an  upper  limit  of 65,000  km 
for  the  loop  size  using  white  light  images  of  the  June  4  flare  show¬ 
ing  loop  footpoints,  assuming  that  such  loops  are  the  same  as 
those  associated  with  the  nuclear  interactions.  The  accelerated 
a/proton  ratio  and  ion  spectral  index  were  next  determined 
simultaneously  using  the  measured  fluence  ratios  of  (1)  the 
a-a  line  complex  and  the  summed  high-FlP  narrow  deexcita¬ 
tion  lines  and  (2)  the  6.129  MeV  160  and  1.634  MeV  20Ne  lines. 
The  a/proton  ratio  was  found  to  be  at  least  0.4  and  the  spectral 
index  was  determined  to  be  3.9  ±  0.2.  Because  the  remaining 
measurable  quantities  depend  to  some  extent  on  the  density- 
height  distribution  of  the  atmosphere,  the  analysis  continued  along 
two  parallel  paths,  each  assuming  one  of  two  possible  atmo¬ 
spheric  models:  the  Avrett  (1981)  model  representing  a  sunspot 
and  a  hybrid  model  combining  the  distribution  above  1 000  km 
derived  from  RUES  SI  X-ray  measurements  (Aschwanden  et  al. 
2002)  with  the  Avrett  model  below,  and  referred  to  as  the 
RHESSU Avrett  model. 

For  the  Avrett  model,  constraints  on  PAS  (1)  and  magnetic 
convergence  (6 )  were  first  attempted  using  the  measured  Doppler 
shift  of  the  4.438  MeV  12C  narrow  deexcitation  line.  But  at  the 
location  of  the  June  4  flare  (74.5°  heliocentric  angle),  the  depen¬ 
dence  of  the  shift  on  the  interacting  ion  angular  distribution  is 
weak  and  significant  constraints  could  not  be  obtained,  other  than 
rejecting  combinations  of  6  and  A,  whose  values  were  simulta¬ 


neously  very  low.  Such  combinations  produce  the  most  redshift 
of  the  line.  Better  constraints  on  6  and  A  were  obtained  by  consid¬ 
ering  the  time  history  of  the  deexcitation  lines.  Several  potential 
time  histories  for  the  accelerated  ion  release  time  history  were  as¬ 
sumed,  and  the  resulting  line  time  histories  were  calculated  for  a 
wide  range  of  combinations  of  6  and  A.  The  calculated  time  his¬ 
tories  were  compared  with  the  measured  time  history  and  \2  was 
used  to  establish  the  quality  of  fit.  Combining  constraints  from 
this  analysis  with  the  weaker  constraints  obtained  from  the  line 
shift  analysis  provided  usable  constraints  on  6  and  A.  We  found 
that  very  low  values  for  both  6  and  A  (corresponding  to  strongly 
downward-directed  interacting  ion  angular  distributions)  and 
very  high  values  for  both  (corresponding  to  very  broad  angular 
distributions)  could  be  rejected.  The  band  of  6-A  combinations 
providing  acceptable  fits  to  both  the  line  shift  and  time  profile 
corresponded  to  mildly  downward-directed  distributions. 

We  next  considered  the  neutron-capture  line.  The  only  remain¬ 
ing  parameter  affecting  this  line  for  which  constraints  had  not  been 
determined  was  photospheric  3He/H.  Neutron-capture  line  time 
histories  calculated  with  varying  values  of  3He/H  and  with  the 
other  parameters  drawn  from  the  acceptable  ranges  previously  de¬ 
termined  were  compared  with  the  measured  time  history.  Accept¬ 
able  fits  to  the  data  were  achieved,  but  the  combined  effect  of  the 
allowed  ranges  for  the  various  parameters  resulted  in  only  an  up¬ 
per  limit  for  3He/H  of  4. 5  x  1 0  5  (3  a),  consistent  with  a  recent 
estimate  of  the  photospheric  3He/H  of  (3.7  ±  0.9)  x  10-5.  We 
then  used  this  recent  3He/H  estimate  to  reverse  the  analysis  and 
place  a  tighter  constraint  on  the  accelerated  ion  power-law  spec¬ 
tral  index  of  3.99  ±  0.06. 

The  final  step  of  the  analysis  was  to  compare  arriving  neutron 
count  rates  calculated  for  the  allowed  ranges  of  parameters  deter¬ 
mined  from  the  previous  analysis  with  neutron  count  rates  mea¬ 
sured  with  OSSE.  We  found  that  the  rates  were  overpredicted 
for  much  of  the  duration  of  the  event,  with  the  largest  excesses 
occurring  early.  Noting  that  early  arriving  neutrons  are  generally 
the  highest-energy  neutrons,  we  calculated  an  arriving  neutron 
count  rate  for  a  power-law  spectrum  having  the  same  spectral 
index  but  with  a  cutoff  above  an  energy  of  125  MeV  nucleon-1, 
which  reduces  the  relative  number  of  high-energy  neutrons.  We 
found  that  an  excellent  fit  to  the  measured  count  rate  was  achieved. 
The  introduction  of  the  cutoff  energy  did  not  significantly  affect 
any  of  the  other  observables  since  they  are  produced  predominantly 
by  accelerated  ions  with  energies  below  125  MeV  nucleon-1. 

Starting  with  the  previously  determined  parameters  that  are  inde¬ 
pendent  of  the  atmospheric  model  (5  =  3 .9  ±  0.2,  and  a/proton  > 
0.4),  we  repeated  the  above  analysis  but  using  the  RHESSU Avrett 
atmospheric  model.  However,  because  of  the  reduced  sensitivity 
of  the  interaction  time  history  at  times  less  than  a  few  seconds  to 
the  interacting  ion  angular  distribution  (see  §  2.3.1),  the  prompt 
deexcitation  line  time  history  analysis  provided  no  useful  con¬ 
straints  on  6  and  A.  We  then  attempted  to  use  the  neutron  capture 
time  history  measurement  (fixing  the  photospheric  3He/H  ra¬ 
tio  at  an  independently  determined  value)  to  constrain  6  and  A  but 
could  only  exclude  6-A  combinations  where  both  parameters  are 
very  small  (corresponding  to  strongly  downward-directed  inter¬ 
acting  ion  angular  distributions).  For  all  acceptable  6-A  combi¬ 
nations,  we  again  found  that  achieving  a  good  fit  to  the  neutron 
arrival  time  history  required  that  the  accelerated  ion  spectrum 
must  cutoff  above  ^125  MeV  nucleon-1. 

We  summarize  what  we  have  learned  about  the  acceleration 
parameters  for  the  1991  June  4  flare.  We  found  that  the  accel¬ 
erated  a/proton  ratio  was  greater  than  0.4.  Assuming  a  power- 
law  form  for  the  spectrum,  we  found  that  the  spectral  index  of  the 
accelerated  ions  from  ~2  to  ~50  MeV  nucleon-1  was  about  4. 
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Above  ^100  MeV  nucleon- 1 ,  the  spectrum  must  steepen  sharply. 
These  parameters  provide  constraints  that  must  be  accounted  for 
by  any  potential  acceleration  mechanism. 

We  have  shown  that  the  measurable  quantities  associated  with 
high-energy  solar  flares  can  depend  in  complex  ways  on  the 
parameters  of  the  transport  and  interaction  model.  Deriving  self- 
consistent  values  for  these  parameters  with  reliable  uncertainties 
can  be  difficult.  Success  requires  measurements  that  cover  a  wide 
range  of  the  observables  and  an  analysis  procedure  that  takes  best 
advantage  of  their  different  parameter  dependences.  We  applied 
such  an  procedure  in  an  analysis  of  the  high-energy  data  available 
for  the  well-measured  1991  June  4  flare.  Even  though  this  flare  is 
one  of  the  best  measured  high-energy  flares,  not  all  of  the  param¬ 
eters  could  be  well  constrained.  Two  factors  contributed  to  this: 
(1)  the  flare  occurred  near  the  solar  limb  so  that  any  line  shift  due 


to  anisotropic  ion  angular  distributions  is  minimal,  and  (2)  be¬ 
cause  of  the  intensity  of  this  flare,  only  data  from  one  off-pointed 
detector  could  be  used  at  the  peak  of  emission,  which  dramatically 
reduced  the  statistical  significance  of  the  time  history  data.  The 
procedure  can  be  applied  to  other  flares  having  a  similar  range  of 
measurable  quantities.  The  successful  self-consistent  prediction  of 
observables  has  been  achieved  within  the  framework  of  our  loop 
transport  and  interaction  model  with  reasonable  values  for  the  as¬ 
sociated  physical  and  acceleration  parameters. 
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APPENDIX 

OSSE  RESPONSE  TO  NEUTRONS 

For  solar  neutron  measurements,  detectors  such  as  OSSE  are  nondiagonal  and  so  do  not  provide  a  good  measurement  of  the  neutron 
energy.  They  do  provide  a  good  measurement  of  the  arrival  time  history  of  the  neutrons.  In  order  to  compare  calculations  of  arriving 
neutrons  with  the  OSSE  measurements,  the  calculated  time-dependent  neutron  spectra  arriving  at  Earth  are  convolved  with  the  OSSE 
neutron  response  to  produce  predicted  neutron  count  spectra.  These  count  spectra  are  then  integrated  over  energy  to  produce  a 
predicted  neutron  arrival  count  rate. 

The  OSSE  instrument  and  its  7-ray  calibration  is  described  by  Johnson  et  al.  (1993).  Here  we  describe  the  neutron  calibration. 
OSSE  consists  of  four  independently  oriented  Nal/Csl  phoswich  scintillation  detectors,  each  surrounded  by  an  anticoincidence  Csl 
shield  and  mounted  in  an  independent  elevation-angle  gimbal,  which  allows  “scanning”  through  192°  of  rotation.  OSSE  can  detect 
both  high-energy  7-rays  and  neutrons  and  the  two  emissions  are  distinguished  by  their  different  pulse  shapes  in  Nal  (Share  et  al. 
1978).  The  identified  energy-loss  events  are  separately  accumulated  into  16  channels  up  to  about  150  MeV  for  7-rays  and  220  MeV 
for  neutrons.  The  separation  is  not  perfect  with  the  result  that  a  fraction  of  neutron  events  are  incorrectly  identified  as  7-rays  and  vice 
versa.  Furthermore,  the  separation  becomes  ineffective  for  neutron  energies  above  ~1  00  MeV  (channel  9). 

Energy-loss  spectra  due  to  all  detected  neutrons  and  to  all  detected  7-rays  can  be  recovered  from  the  separated  spectra  by  using  the 
separation  constants  (i.e.,  the  energy-dependent  fractions  of  neutron  events  that  are  identified  as  neutrons  and  as  7-rays,  and 
analogously  for  the  7-ray  events)  developed  by  Del  Signore  (1995).  These  constants  were  derived  for  23-100  MeV  neutron  energy 
losses  (channels  3-8)  by  assuming  that  the  >16  MeV  emission  at  the  peak  of  the  1991  June  6  flare  (  which  did  not  produce  neutrons) 
was  purely  7-rays,  while  the  emission  near  the  end  of  the  1991  June  4  event  (after  all  evidence  for  7-ray  emission  has  subsided)  was 
purely  neutrons.  We  do  not  consider  neutron  energy  losses  <23  MeV  (channels  1  and  2)  since  at  these  energies  the  neutron  response  is 
falling  rapidly  and  the  count  rates  are  dominated  by  the  lower  energy  7-rays  arriving  simultaneously  in  the  June  4  flare.  Above  100  MeV 
(channel  9),  where  the  total  neutron  and  7-ray  spectra  are  not  recoverable  from  the  separated  spectra,  we  find  for  the  June  4  flare  that 
whether  these  events  are  assumed  to  be  all  neutrons  or  all  7-rays,  the  resulting  total  neutron  count  rate  is  affected  only  early  in  the  event 
(when  such  high-energy  neutrons  are  arriving)  and  then  only  insignificantly.  We  use  the  separation  constants  determined  at  100  MeV  for  all 
>100  MeV  energy  losses. 

A  prototype  OSSE  detector  was  exposed  to  neutrons  at  the  Indiana  Cyclotron  (Jenkins  et  al.  1991)  and  an  instrument  model, 
consisting  of  energy-deposited  spectra  for  several  incident  neutron  energies,  was  developed  (T.  Jenkins  2000,  private  communication). 
When  these  energy-loss  spectra  are  integrated,  the  resulting  total  OSSE  neutron  response  showed  a  strong  turnover  at  neutron  energies 
above  80  MeV.  While  a  turnover  is  expected  (due  to  self-veto  of  high-energy  neutrons  depositing  energy  in  both  detectors  and 
shields),  it  is  not  expected  to  be  so  severe  and  the  turnover  onset  energy  should  be  above  ~  1 00  MeV,  as  derived  by  Cooper  et  al.  (1 985) 
for  the  similar  SMM GRS  Nal/Csl  detectors.  In  addition,  the  total  OSSE  neutron  effective  area  at  these  energies  appears  to  have  been 
underestimated  since  it  is  less  than  that  of  the  smaller  SMM  GRS  Nal/Csl  detectors.  To  correct  the  response,  we  used  the  Nal  neutron 
response  measurements  (20  MeV  threshold)  of  Stanislaus  et  al.  (200 1 )  to  renormalize  each  neutron  energy-loss  spectrum,  maintaining 
their  relative  energy  dependences  as  derived  by  Jenkins.  Above  the  neutron  energy  of  1 50  MeV,  we  introduced  a  roll-off  of  the  Stanislaus 
effective  areas  (which  continue  to  rise  with  neutron  energy),  using  a  neutron-energy  dependence  similar  to  that  of  the  SMM  GRS  given  by 
Cooper  et  al.  (1985).  The  change  in  the  neutron  count  rate  predicted  by  this  modified  response  compared  to  the  count  rate  predicted  by  the 
original  response  established  by  Jenkins  is  small  since  neutrons  >100  MeV  do  not  contribute  significantly  to  the  total  detected  count  rate 
and  is  within  the  estimated  30%  overall  uncertainty  of  the  neutron  response  that  we  have  included  in  the  June  4  analysis. 
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